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The effects of oscillating flow on the pressure force normal to the 
chord of a symmetrical airfoil were investigated experimentally employ- 
ing a remote pressure transducer to measure the instantaneous pressure 
distribution . 

An open circuit wind tunnel having a set of rotating shutter blades 
located down stream of the test section was used to produce the oscil- 
lating flow. Electrical signals analogous to the free stream velocity and 
surface pressure were recorded simultaneously on separate tracks of a 
magnetic tape. The recorded data were converted to digital representation, 
and numerical techniques utilized to evaluate the spectral composition of 
the measured pressure distribution, from which the normal force was 
ca leu lated . 

It was found that the magnitude of the total normal force at high 
angles of attack is significantly greater in oscillating flow than in 
steady flow and is frequency dependent ; while at low angles of attack 
no significant differences were observed. Moreover it was found that 
higher order harmonics of the fundamental free-stream frequency constitute 
a significant fraction of the normal force, and these fractions are also 
frequency dependent. The observed results are not adequately predicted 
by quasi-steady aerodynamic analysis. 
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I . INTRODUCTION 



The aerodynamic forces associated with unsteady flows have attracted 
great interest in recent times. With the rapid development of the heli- 
copter and advances in gas turbine technology, more detailed knowledge 
of unsteady phenomona is required than can be obtained from quasi-steady 
flow ana lys is . 

Experimental investigations of unsteady flow problems have been 
hampered by lack of adequate facilities, i.e. lack of wind tunnels that 
are capable of establishing a controlled unsteady flow. Consequently 
the majority of previous investigations have been concerned with the 
starting process or have utilized programmed model motion rather than 
an oscillating flow [Ref. 1, 2, and 3]. Analytical investigations of 
the problem have been restricted to linearized solutions that are not 
applicable for airfoils at large angles of attack [Ref. 4 and 5]. 

The objective of the work reported here was to measure the time- 
dependent pressure distribution on an airfoil at constant angle of 
attack in a controlled oscillating stream and from these measurements 
to calculate the time-dependent forces acting on the airfoil surface. 
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II. BASIC APPROACH 



A. ALTERNATIVES AVAILABLE FOR MEASUREMENT 

There are two obvious methods one might use to measure the instan- 
taneous pressure distribution on the surface of an airfoil model in an 
unsteady flow. The most direct technique would be to mount many trans- 
ducers in the model in such a manner that the pressure at discrete 
points on the surface could be measured directly. The difficulties 
involved in applying this method are the cost of the miniature trans- 
ducers required and the possibility that the dynamics of the test 
installation would influence the response of the transducers. 

A less direct method would be to connect remote transducers to 
pressure taps on the surface of the airfoil with tubing. This would 
eliminate the dynamic interference likely to be encountered with the 
more direct method, but would add to the complexity of the system 
because of the transfer functions associated with the tubing. 

In the case of oscillating flows, however, both the free stream 
velocity and the pressure at a point on the surface of the airfoil are 
periodic functions of time having some fixed phase relation between 
them. Knowledge of the phase relations existing between the free 
stream velocity and the pressure at various points on the surface of 
the airfoil, assuming that this phase relation remains constant, would 
allow the pressure distribution to be reconstructed from temporal 
records of the pressures and velocities taken over different intervals 
of time. A remote transducer then could be utilized to measure pres- 
sure and, if the leads from each tap were identical, only a single 
transfer function would have to be known. 
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B. BRIEF DESCRIPTION OF EXPERIMENTAL METHOD 



The technique utilized in this investigation was based on the last 
alternative described above. Electrical signals analogous to the free 
stream velocity and surface pressure at a point were recorded simul- 
taneously on separate tracks of magnetic tape. The pressure analogues 
recorded were the electrical output of a single remote transducer con- 
nected in turn by tubing to each of the pressure taps on the surface of 
the airfoil. The velocity analogues recorded were the output of a 
linearized hot-wire anemometer. 

The transfer function of the tubing was experimentally measured and 
applied in digital form to the Fourier Transform of the digitized 
pressure analogue in order to obtain the Fourier Transform of the 
pressure at the surface of the airfoil. The inverse of the pressure 
Fourier Transform was then taken to obtain the pressure- time relation. 
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III. EXPERIMENTAL EQUIPMENT 



A. WIND TUNNEL 

1 . Genera 1 Description 

The experimental work was conducted in the low-speed, oscillating 
flow wind tunnel located in the Aeronautics Laboratories of the Naval 
Postgraduate School, This wind tunnel is of open circuit design, with 
a 24-inch square by 223-inch long test section. A plan view of the 
tunnel is presented in Figure 1. The tunnel inlet is eight feet square, 
resulting in a 16:1 contraction ratio. Three high solidity screens 
located in the inlet section just upstream of the nozzle produce measured 
free stream turbulence intensities of 0.261 to 0.413 percent for the 
velocities encountered in the present work. 

The wind tunnel drive consists of two Joy Axivane Fans in series, 
each of which has an internal, 100 horsepower, direct connected, 1750 rpm 
motor. The fan blades are internally adjustable through a pitch range of 
25 to 55 degrees, providing a wide operating base. Two sets of variable 
inlet vanes, located immediately upstream of each fan, are externally 
operated to provide control of test section velocity. These vanes are 
of multileaf design, and preswirl the air in the direction of fan 
rotation to reduce fan capacity. The total range of tunnel velocity is 
from 10 to 250 feet per second. 

2 . Rotating Shutter Valve 

Two fundamental methods of creating an oscillating flow environ- 
ment have been employed in the past. Both Nickerson [Ref. 6] and Hori 
[Ref. 7] introduced oscillations by oscillating their models in a steady 
flow environment. This method severely restricts the range of attainable 
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PLAN VIEW OF WIND TUNNEL 



frequencies because of mechanical complications, and also introduces 
measurement difficulties. The other approach is to actually oscillate 
the flow over a stationary model. Hill [Ref. 8] used a sliding shutter 
to impose oscillations on the free stream but was restricted by mechanical 
limitations to low frequencies. 

The most successful method of obtaining an oscillating flow with 
large ranges of frequency and amplitude was that employed by Karlsson 
[Ref. 9], and later by Miller [Ref. 10] in his investigation of trans- 
ition. A rotating shutter valve, immediately downstream of the test 
section, is employed to superimpose a periodic variation of velocity on 
the mean flow. The method used in the present investigation is virtually 
identical to that employed by Miller. The shutter valve consists of four 
horizontal steel shafts equally spaced across the test section. The 
shafts are slotted to accommodate flat blades of various widths, forming 
a set of four butterfly valves spanning the test section. Figure 2 is a 
photograph of the shutter valve. Each blade is driven from its immediate 
neighbor by means of a timing belt and pulley arrangement. The bottom 
shaft is driven by a five horsepower variable speed electric motor, 
through a timing belt and pulley. An intermediate shaft between the 
motor and shutter valve permits a wide variety of pulley ratios. This 
drive arrangement provides a frequency range of two cycles per second to 
the first critical frequency of 933 cycles per second. The electric 
motor presently in use, however, restricts the oscillation frequency to 
a maximum of 240 cycles per second. The amplitude of oscillation is con- 
trolled by blade width. Test section closure may be varied from 25 to 100 
per cent. The resulting amplitude of oscillation of test section velocity 
is a function of frequency, mean velocity and pressure gradient. In this 
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GjRE l photograph of rotating shutter valve 





investigation, blades producing 50.0 per cent closure were used, resulting 
in an amplitude range of from' 3 to 40 per cent of the local mean free 
stream velocity. 

3 . Test Section 

The wind tunnel test section is shown in Figure 3. Continuous 
pieces of two-inch thick aluminum, 24 inches wide and 223 inches long, 
form the upper and lower test section walls. Each of the side walls 
consists of three two-inch thick panels of stress -re lieved Lucite. For 
this investigation, the central side wall panel on the opposite side of 
the tunnel from the control console was replaced with two-inch thick 
plywood to facilitate the mounting of instrumentation. The Lucite panels 
on the console side of the test section are hinged and may be raised 
hydraulically, providing access to the test section. The heavy con- 
struction of the test section is intended to minimize deflections 
induced by rapid changes in static pressure. 

Figure 4 is a typical test section velocity profile. Velocity 
variation is less than one per cent from the mean to within three inches 
of any wall. 

Figure 5 is an overall photographic view of the wind tunnel. 

B. MODEL 

The airfoil model used in this investigation was that employed by 
Allen [Ref. 11]. It is a NACA 63-010 section modified by straight line 
fairing from 60 per cent chord to the trailing edge. The model has a 
24 inch span and a constant 6 inch chord. Twenty three pressure taps 
are located chordwise across the upper surface and two on the lower 
surface at midspan. Figure 6 is a schematic drawing of the model showing 
the location of the pressure taps. 
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FIGURE 6. PRESSURE TAPS ON WIND TUNNEL MODEL 



The model was mounted spanwise across the test section on a mechanism 
which allowed measuring the angle of attack about the midchord line. The 
pressure tap leads were brought through the tunnel walls with stainless 
steel tubing. Figure 7 shows the model mounted in the tunnel test section. 

C. PRESSURE MEASURING SYSTEM 
1 . Pressure Transducers 

Limitations on size and expense lead to the selection of a remote 
pressure transducer designed by Professor L. V. Schmidt of the Naval 
Postgraduate School. The design is shown in Figure 8. Two transducers 
were built, one to measure pressures on the model and the other to be 
used as a reference. 

The transducer system consisted basically of a Bently Detector 
System, a 0.003 inch thick annealed brass diaphragm mounted in an 
aluminum housing, and plastic and steel tubing enabling the pressure on 
the airfoil to be transmitted to the transducer diaphragm cavity. The 
diaphragm thickness was selected so as to provide a useable instrument 
in the pressure range of + 0.75 psia. The disassembled transducer is 
shown in Figure 9. 

The Bently Detector System provides a voltage signal which is 
linear with respect to the distance change detected by its probe. For 
small pressure differentials across the brass diaphragm, the deflection 
is proportional to pressure. By putting a known pressure differential 
across the diaphragm, the transducer can be calibrated yielding a 
linear calibration curve for small differentials. The details of this 
procedure are given in Appendix A. 

The stainless steel and plastic tubing act as transmission lines 
from the pressure taps to the transducer. The unsteady pressure components 
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FIGURE 8. SECTIONAL DRAWING OF PRESSURE TRANSDUCER 
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FIGURE 9 PHOTOGRAPH OF TRANSDUCER COMPONENTS 




undergo frequency-dependent phase shift and attenuation through the 
tubing. A dynamic calibration (Appendix A) was performed at 5 hertz 
increments and intermediate values at 1 hertz increments were obtained 
by linear interpolation and extrapolation. The results may be seen in 
Figures 10 and 11. 

2 . Power Supplies and Amplifier 

A Power Design Transistorized Power Supply, Model 205 1R was used 
to power the Bent ly Proximeter. The Proximeter output at zero pressure 
differential was in the neighborhood of + 5.73 Volts as installed in the 
pressure measuring system. A Kepco Model OPS 21-1 Operational Power 
Supply/Amplif ier was used to amplify the Proximeter signal. The output 
of the operational amplifier was adjusted to zero volts at zero pressure 
differential by means of a suitable bias voltage. The gain of the 
operational amplifier at a given feedback resistance was verified to be 
constant with input frequency over the frequency range of interest. The 
phase shift due to the operational amplifier was a constant -180 degrees 
over the same frequency range. Figure 12 is a schematic diagram of the 
pressure measuring system and Figure 13 is a typical static calibration. 
Figure 14 shows the pressure measuring system installed. 

D. HOT WIRE ANAMOMETER 

A constant temperature, transistorized, hot wire anamometer was used 
to measure the free stream velocity. The electrical output of the hot 
wire anamometer was proportional to the velocity. A micromanometer 
connected to a pitot-static tube was used to measure the mean free 
stream velocity and to statically calibrate the hot wire anamometer. 
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FIGURE 10. TRANSDUCER GAIN AS A FUNCTION OF FREQUENCY 



Frequency, Hertz 
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FIGURE II. TRANSDUCER PHASE ANGLE, 
AS A FUNCTION OF FREQUENCY 
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FIGURE 12. SCHEMATIC DIAGRAM OF PRESSURE MEASURING SYSTEM. 





FIGURE 13. TYPICAL STATIC CALIBRATION 
CURVE OF PRESSURE MEASURING SYSTEM 
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E . TAPE RECORDER 



The outputs of the hot wire anamometer and pressure transducer were 
recorded on an Ampex CP*- 100 Portable Instrumentation Magnetic Tape 
Recorder/Reproducer. Typical record and reproduce calibration values 
for the recorder at a tape speed of 3-3/4 inches per second are shown 
in Table I. 

TABLE I 

TYPICAL TAPE RECORDER CALIBRATION 
Record 

D.C. INPUT (VOLTS) FREQUENCY (HERTZ) 

1.414 9448 

0 6748 

-1.414 4047 



Reproduce 



FREQUENCY INPUT (HERTZ) D.C. OUTPUT (VOLTS) 

9450 1.420 

6750 0.0 

4050 -1.409 

F. MISCELLANEOUS 

The velocity and pressure analogues were monitored on a Tektronix Dual 
Beam Oscilloscope. Pictures of representative velocity and pressure 
outputs were recorded photographically from oscilliscope traces. 

The shutter valve frequency was measured with a magnetic pickup 
located outboard of the upper shutter valve blade shaft. The output of 
the pickup was read on a decade counter. A schematic diagram of the 
complete data aquisition system is shown in Figure 15. 
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FIGURE 15. SCHEMATIC DIAGRAM OF DATA 
ACQUISITION SYSTEM ' 
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IV. PROCEDURES 



Prior to each run, a static calibration of the transducer was made in 
order to ensure that no drift had occurred since the previous calibration. 

The desired free stream mean velocity was established in steady flow 
and the hot wire anamometer set to give the desired output at the free 
stream mean velocity. The oscillations were then adjusted to the required 
frequency . 

The angle of attack vernier was calibrated by comparing the trans- 
ducer output of the taps on the lower surface with those at corresponding 
locations on the upper surface. The angle of attack vernier was accurate 
to + .05 degrees when calibrated on this basis at zero degrees angle of 
attack. 

Apparent regularity of the free stream wave-form on the oscilloscope 
was used to establish the desired oscillating frequencies. Subsequent 
Fourier analysis of the free stream velocity wave-form showed this pro- 
cedure to be adequate in establishing clean wave-forms. 

The airfoil was placed at the desired positive angle of attack and 
recording began. The output of the transducer connected to each pressure 
tap in turn was recorded on tape simultaneously with the output of the 
hot wire anamometer. Since the airfoil was symmetrical, the pressure on 
the lower surface was assumed to correspond to the pressure on the upper 
tapped surface at negative angle of attack. The airfoil was rotated to 
a negative angle of attack in order to obtain the pressure distribution 
on the lower surface. 
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V. DATA ANALYSIS 



A. ANALOGUE TO DIGITAL CONVERSION 

The recorded data were converted from analogue to digital repre- 
sentation on the Naval Postgraduate School’s Hybrid Computer. Details 
of the data conversion process are given in Appendix B. 

The analogue signals were filtered through two low-pass R. C. filters 
in series prior to digitization. The filter was used to improve signal 
to noise ratios of the analogue data. Preliminary investigations had 
shown that some high frequency noise was picked up in the data recording 
process from both the transducer and the tape recorder. ( M High frequency*', 
for the purposes of this investigation, means frequencies greater than 500 
hertz.) The filters were built around the existing high gain amplifiers 
on the analogue side of the hybrid. The filter transfer function was 
measured at 10 hertz increments and intermediate values at 1 hertz 
increments obtained by linear interpolation. Figures 16 and 17 give 
the transfer function of the filtering system. 

All sampling was done at the rate of 1024 samples per second. The 
analogous pressure and velocity signals associated with each pressure tap 
were sampled for 4 seconds for a total sample size of 4096 samples per 
digital record. 

Prior to each run a D. C. reference voltage measured with a digital 
voltmeter was recorded. These reference voltages were then digitized as 
a check on the tape recorder and digitizing process. Typical values of 
voltmeter measured D. C. voltages and the corresponding digital mean 
value are shown in Table II. 
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FIGURE 17. FILTER PHASE ANGLE AS A FUNCTION OF FREQUENCY 



TABLE II 



RECORDED D.C. VOLTAGE AND CORRESPONDING DIGITAL MEAN VALUE 



RECORDED (VOLTS) 



DIGITAL MEAN VALUE (VOLTS) 



+ 1.414 



+1.409 

0.995 

0.496 

0.104 

0.002 

- 0.100 

-0.496 

-0.994 

-1.415 



1.0 

0.5 

0.1 

0.0 

- 0.1 

-0.5 

- 1.0 



-1.414 



B. TAPE CONVERSION 

The arrays of digitized data obtained from the analogue to digital 
conversion process were stored on 7 track, 200 bit per inch magnetic 
tape in 24 bit words. A conversion to 9 track, 800 bit per inch tape 
and 32 bit words was necessary in order to analyze the data on the IBM 
360 Series Computer. The computer program used to make the conversion 
is included as an appendage to this thesis under the title "TAPE CON- 
VERSION PROGRAM". 

C. ADDITIONAL ASSUMPTIONS 

Ideally, phase information between the free stream velocity and 
pressure analogues at a given tap would have been preserved during the 
analogue to digital conversion process. Because of equipment limitations 
discussed in Appendix B, this was not possible. However, measurements 
from photographs of oscilloscope traces of the velocity and pressure 
analogues corresponding to given pressure taps indicated that the 
phasing between the free stream velocity and the pressure on the surface 
of the airfoil was essentially independent of the surface coordinate. 

This measured time difference between a negatively-sloped zero crossing 
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of the corresponding pressure analogue was therefore taken as constant 
for a given run. 



D. NON-DIMENS IONALIZATION 

1. Pressure Coefficients 

The pressure data corresponded to the difference between atmos- 
pheric pressure and the surface pressure. It was therefore convenient 
to define a non-dimensional pressure coefficient referenced to atmos- 
pheric pressure as: 

P(x,y,t)-P 

C p (x,y,t) = z 

q 

where : 



P(x,y,t) = pressure at time t at a point (x,y) on the surface of 
the airfoil, 

P = atmospheric pressure, and 

A 

q = mean free stream dynamic pressure. 

The mean pressure coefficient, Cp(x,y), is then defined as: 

P(x,y)-P 

C p (x,y) = and 

q 

the "unsteady pressure coefficient" is defined as: 



C (x,y,t) 



p( x >y> t ) 
> 

q 



where : 



p(x,y,t) = "unsteady pressure" at time t at a point (x,y) on the 
surface of the airfoil. 

2. Normal Force Coefficients 



The time dependent force per unit span acting normal to the chord 
of a two-dimensional airfoil is approximated by: 
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c 



N(t) =J [ P(x , -y , t )-P(x ,y , t ) ] dx 
0 

if the airfoil is symmetrical. The normal force coefficient can be 

defined in terms of the non-dimensional pressure coefficient as: 

1 

c N (t) = J [C p (x,-y,t)-C p (x J y,t)) d§ 

0 



where : 

§ = x/c. 

In terms of the mean and unsteady pressure coefficients this becomes: 

1 1 

C N (t) = / [C p (S f Ti) ' V §,Tl)1 dZ + I dc 

0 0 



which leads to: 






+ MU 

qc qc 




+ c n (t) 



where : 



C x , = mean normal force coefficient 
N 

c^(t) = ’’unsteady normal force coefficient". 

3 . Center of Pressure 

The center of pressure may be defined in terms of non-dimensional 

parameters in a manner similar to that used above. The result is: 

1 

j [C p (S,-Tl,t) - C p (S,r,t)] § d§ 

F = 5 

C -P- C N (t) 

E. COMPUTATIONAL PROCEDURES 

The numerical techniques used in the calculations are discussed in 
Appendix C. A flow chart of the computer program used to analyze the 
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digitized data is shown in Fig. 18. A listing of the program is included 
as an appendage to this thesis. 
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FIGURE 18. ABBREVIATED COMPUTER PROGRAM FLOW CHART 
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VI. RESULTS AND DISCUSSION 



A. CONDITIONS INVESTIGATED 

Table III lists the freestream conditions and angles of attack at 
which data was taken. A "Run" consisted of recording the pressure at 
each tap for 30 seconds with the airfoil at a given angle of attack in 
a free stream velocity with a given mean velocity, amplitude of oscil- 
lation, and frequency of oscillation. 

The steady flow Runs were made and the pressure distribution and 
normal force coefficients presented in order to provide a basis for 
comparison with the non-steady results. Power spectral densities of 
the turbulent pressure fluctuations in steady flow were calculated to 
investigate the possibility of any harmonic coupling that might occur 
between the forced oscillations of the free stream and the vortex 
shedding phenomona observed in steady flow. 

B. STEADY FLOW RESULTS 

The measured pressure distributions on the airfoil at angles of 
attack of 5, 10, and 15 degrees in a steady flow of approximately 100 
feet per second are shown in Fig. 19 through 21. Figure 22 shows the 
resulting normal force coefficients. The lift curve of a NACA 63-009 
airfoil section is included for comparison [Ref. 12]. 

Power spectral densities of the turbulent pressure fluctuations at 
representative taps are shown in Fig. 23 through 32. These power 
spectrums and subsequent power spectrums shown have been normalized with 
respect to their mean squares. 

At angles of attack of 0 and 5 degrees a distinct periodic component 
in the neighborhood of 2 1 to 22 hertz was discernable at all chordwise 
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locations investigated. The corresponding Strouhal Number taken with 
respect to chord and the frequency of the periodic component is : 

N ST = °* 105 - 

At 10 degrees angle of attack the spectrum of the fluctuating 
pressure on the upper surface appears to be broadened over the lower 
frequencies. Periodic components are discernable along the lower sur- 
face at about the same frequency as was noted at lower angles of attack. 

At 15 degrees angle of attack the airfoil was fully stalled as 
evidenced by the pressure distribution. Figures 30 through 32 show a 
distinct periodic component between 21 and 22 hertz just as the runs at 
lower angle of attack. This periodicity was visible at all pressure 
taps on the surface. The tendency of the spectrum to spread more 
towards lower frequencies at downstream stations on the upper surface 
is apparent when Fig. 30 and 31 are compared. 

Figures 33 and 34 show the mean square of the fluctuating pressure 
coefficients for Runs 3 and 4 plotted as a function of chord. Discussion 
of the possible significance of the resulting distribution will be 
deferred to the unsteady flow section. 

Noise at 60 hertz was filtered from the spectral and mean square 
distributions shown in Fig. 23 through 34 by setting the computed value 
of the spectrum at 60 hertz equal to the average of tbe values at 59 and 
61 hertz. A typical unfiltered spectrum is shown in Fig. 35. 

C. UNSTEADY FLOW RESULTS 

1 . Run 5 

Run 5 was made at a fundamental frequency of 94 hertz, an angle 
of attack of 15 degrees, and a free stream velocity with mean of approxi- 
mately 100 feet per second and amplitude ratio of 0.085. Figures 36 and 
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37 show typical oscilloscope traces of the velocity and pressure analogue 
wave-forms at selected pressure taps for this Run. 

It is apparent from these photographs that the time between a 
negatively sloped zero crossing on the velocity analogue wave and a 
positive sloped zero crossing on the pressure analogue wave is essentially 
constant at a given tap even though the wave-form is obviously not the 
same for every period. These same pictures show that the change in the 
above time-difference with chordwise location is small. 

The time-difference mentioned above was measured at two sample 
time-increments from photographs similar to those shown in Fig. 36 and 
37, but on an expanded time scale. Expressed in terms of a phase angle 
at the fundamental frequency of 94 hertz, the two sample time-increments 
are equivalent to 66 degrees. This measurement served to establish the 
phasing that was used in reducing the data. 

Figure 38 shows the power spectrum of the oscillating free 
stream velocity for Run 5. The spectral quality of the oscillations 
produced by the tunnel is immediately apparent from the manner in which 
the fundamental frequency dominates the spectrum. 

Figures 39 through 64 are power spectral densities of the unsteady 
pressures associated with the indicated chordwise locations. They were 
selected to indicate the complete trend of the pressure spectrum on the 
surface of the airfoil. The spectrums appear to be essentially the same 
along the upper surface to about 30 per cent chord at which point the 
spectral component at the fundamental harmonic begins to decrease and the 
second and third order harmonic components increase. These effects appear 
to be maximized at or near the mid-chord point where the second order 
harmonic dominates the unsteady pressure. (It is unfortunate that the 
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spacing between pressure taps near the mid-chord point increases but the 
observed pressure spectrums were not anticipated before the model was 
constructed.) Less than 20 per cent chord downstream from the mid-chord 
point the spectrum again takes on the appearance noted near the leading 
edge . 

Figure 65 shows the measured mean square of the unsteady pressure 
coefficient at each tap. There is a decrease in the mean square as mid- 
chord is approached and a definite increase as the downstream edge is 
approached, implying some sort of spatial correlation between the mean 
square of the pressure fluctuation and chordwise distance downstream. 

Figure 66 shows the calculated time dependent normal force 
coefficient and the free stream velocity plotted with common time origins 
as determined by the previously described photographic measurement. 

Figure 67 is the power spectrum of the calculated unsteady normal 

force . 

Figure 68 shows the mean pressure distribution as calculated for 
Run 5. Figures 69 and 70 are photographs of the time dependent pressure 
distribution displayed on an IBM 2250 Visual Display Unit. The computer 
program used to create the display is attached as an appendage under the 
title "DIS PLAY PROGRAM' 1 . 

Figure 7 1 is the power spectrum of the calculated unsteady 
center of pressure for Run 5. 

2 . Run 7 

Run 7 was made at a fundamental frequency of 11 hertz, an angle 
of attack of 15 degrees, and a free stream velocity with a mean of approxi- 
mately 100 feet per second and amplitude ratio of 0.072. Figures 72 and 
73 show typical oscilloscope traces of the velocity and pressure analogue 
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wave-forms at selected pressure taps for Run 7. Note the regularity of 
the pressure analogue wave-forms with time for a given tap and the 
similarities of the wave-forms for the different taps. The time dif- 
ference between a negatively sloped zero crossing on the velocity 
analogue wave-form and a positively sloped zero crossing on the pressure 
analogue wave-form was measured at approximately 24 sample time -increments . 
At the fundamental frequency of 11 hertz, 24 sample time-increments is 
approximately given in terms of a phase angle as 93 degrees. This measure- 
ment was used as described previously to establish the phasing between the 
free stream velocity and pressure in data reduction. 

The power spectral density of the unsteady free stream velocity 
for Run 7 is shown in Fig. 74. The quality of the oscillation should 
again be noted. 

Power spectrums of the unsteady pressure associated with repre- 
sentative taps are shown in Fig. 73 through 80. The wave-form similarity 
noted in the oscilloscope pictures is also apparent in the power spectrums. 
The similarity of the spectrums to each other should also be noted. 

Figure 81 gives the mean square of the unsteady pressure coefficient 
associated with each tap as a function of the non-dimensional chord. 

The unsteady velocity and normal force coefficient are given as 
functions of time from a common time origin in Fig. 82 and 83. The power 
spectrum of the unsteady normal force is given in Fig. 84 and the power 
spectrum of the center of pressure in Fig. 85. Figure 86 shows the mean 
pressure distribution on the airfoil for Run 7. 

D. COMPARISON OF RESULTS AT DIFFERENT FREQUENCIES 
1 . Normal Force and Center of Pressure 

The most obvious difference between the results presented is in 
the mean normal force coefficients. This is illustrated in Fig. 87 which 
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shows the mean normal force coefficient at 15 degrees angle of attack as 
a function of frequency. Recall that the mean velocities for the steady 
flow Run and the oscillating Runs 5 and 7 were all about 100 feet per 
second. Run 5 had a fundamental frequency of 94 hertz while Run 7 had 
a fundamental frequency of 11 hertz. The amplitude ratio for the two 
oscillatory Runs was about the same, that for Run 5 being 0.085 and that 
for Run 7 being 0.072. 

The harmonic composition of the unsteady normal force coefficient 
for the two oscillatory Runs under consideration was also much different. 

A comparison of Fig. 67 and Fig. 84 shows that the second harmonic of 94 
hertz is much more predominant than the second harmonic of 11 hertz. 

Center of pressure calculations also reflect the difference in 
harmonic composition that was found in the normal force results. The 
unsteady center of pressure at 94 hertz is dominated by the second harmonic 
and the first and third harmonics are about the same, while at 11 hertz 
little difference between the spectrum of the unsteady normal force and 
unsteady center of pressure was noted. 

2 . Pressure Distributions 

The spectral composition of the unsteady pressure reflects the 
difference noted in the normal force and center of pressure results. The 
presence of higher order harmonics is apparently not just related to 
increasing frequency. Allen [Ref. 11] showed oscilloscope traces of the 
pressure analogues associated with the same airfoil used in this investi- 
gation at a fundamental frequency of 128 hertz and about the same mean 
and mean square velocity. His traces did not give an indication of any 
significant higher order harmonics. 

There are apparent spatial periodicities evident in the instan- 
taneous pressure distributions shown in Fig. 69 and 70 for Run 5. The 
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mean pressure distribution for Run 7 shown in Fig. 86 also clearly 
indicates a spatial periodicity. The mean square distributions of the 
unsteady pressure plotted in Fig. 65 and 81 show a similar trend. 

A hint as to the cause of these periodic trends may be obtained 
by considering the chordwise distribution of the mean square of the 
pressure fluctuation in steady flow. One would normally expect vortices 
shed at or near the separation point on the airfoil to immediately begin 
to dissipate as they travel downstream from their point of origin. A 
resulting decrease in the mean square of the fluctuations in velocity 
and pressure would be expected with chordwise distance downstream as is 
shown in Fig. 33 for the airfoil at 10 degrees angle of attack in steady 
flow. The mean square of the pressure fluctuation on the airfoil at 15 
degrees angle of attack in steady flow shows a definite increase toward 
the trailing edge as shown in Fig. 34. This suggests that the spatial 
periodicity observed in the unsteady runs was a result of the combination 
of the oscillatory effects and the effects of the stalled airfoil. 

The instantaneous pressure distributions shown for Run 5 in Fig. 

69 and 70 also suggest that the second harmonic component of the normal 
force could be due to a delay of separation near the leading edge associ- 
ated with the minimum of the velocity or change in sign of the acceleration 
of the free stream. Reference to the instantaneous pressure distribution 
associated with the minimum of the velocity demonstrates the more favorable 
spatial pressure distribution that occurs at that point on the velocity 
wave form. 

E. COMPARISON OF UNSTEADY FLOW RESULTS WITH STEADY FLOW ANALYSIS 

Steady flow results have traditionally been used to estimate aero- 
dynamic forces acting on airfoils in unsteady environments. Such analysis 
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usually involves assuming that the lift coefficient of an airfoil section 
is constant at a given angle of attack within a given range of Reynolds 
Number . 

The normal force coefficient found for the airfoil used in this 
investigation at 15 degrees angle of attack and a Reynolds Number based 
on chord of approximately 3 x 10^ was 0.681. Assuming the normal force 
coefficient to be constant for a given angle of attack leads to the fol- 
lowing estimates for the mean normal force coefficient and mean square of 
the unsteady normal force coefficient for the same airfoil at 15 degrees 
angle of attack in an oscillating flow with amplitude ratio e = 0.08: 

~ 0.683 
s 

%’ C „ (2e 2 +0=0.0128 

s 

where C„ is the normal force coefficient in steady flow. 

N 

s 

Note that assuming a constant normal force coefficient leads to 
results that depend only on the amplitude ratio of the oscillating 
free stream and are independent of frequency. 

The experimentally measured mean normal force coefficients for the 
airfoil at 15 degrees angle of attack in a free stream oscillating at 11 
and 94 hertz were 0.851 and 1.401 respectively. The mean squares of the 
unsteady normal force coefficients for the two cases were 0.075 and 0.100 
respectively. The difference between the experimentally derived values 
and the values predicted by assuming a constant normal force coefficient 
at a given angle of attack are considered by the author to be significant. 

F. SUMMARY OF RESULTS 

Figure 87 gives the measured mean normal force coefficient as a 

function of frequency for the airfoil at an angle of attack of 15 degrees 
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in a free stream with a mean velocity of approximately 100 feet per 
second and amplitude ratio of approximately 0.08. Figure 88 gives the 
measured mean center of pressure under the same conditions. 

Due to lack of computing time and a machine malfunction in the 
digitization process, the data taken for Runs 6 and 8 were not reduced. 
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VII. CONCLUSIONS 



From the results presented, the following conclusions may be drawn 
for the conditions investigated: 

1. The magnitude of the pressure force acting normal to the 
chord of an airfoil at high angle of attack in oscillating flow is 
frequency dependent and is significantly greater than that observed in 
steady flow and is not adequately predicted by quasi-steady analysis. 

2. Higher order harmonic components of the fundamental free 
stream frequency of oscillation are of significant order of magnitude in 
the normal force and their relative magnitude is frequency dependent. 

3. Distinct periodic components are observable in the pressure 
acting on an airfoil in steady flow. 
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TABLE III 



EXPERIMENTAL CONDITIONS FOR WHICH DATA WERE TAKEN 



Run No. 


Angle 

of 

Attack 

(degrees) 


Mean 
Ve locity 
(ft/ sec) 


Frequency 

of 

Oscillation 
(hertz ) 


Amplitude 

Ratio 

of 

Oscillation 


i 


0 


100 


steady 


-- 


2 


5 


100 


steady 


-- 


3 


10 


100 


s teady 


-- 


A 


15 


100 


steady 


-- 


5 


15 


100 


94 


0.08 


6 


15 


150 


94 


0.08 


7 


15 


100 


11 


0.08 


8 


15 


70 


9 


0.25 
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o NACA 63-009 STEADY FLOW LIFT CURVE Nr E = 3x I0 6 Cref. ] 
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MEAN PRESSURE DISTRIBUTION ON AIRFOIL 
> =10 degrees STEADY FLOW 
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MEAN PRESSURE ON AIRFOIL 
CX = 15° STEADY FLOW 
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£,= 0.0 MEAN SQUARE 
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FIGURE 25. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 2 
£»0.0 MEAN SQUARE =0.00334 
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Frequency , Hertz 

FIGURE 26. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 2 
£ = O.I5 UPPER SURFACE MEAN SQUARE =0.00993 
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FIGURE 27. PSD OF UNSTEADY PRESSURE COEFFICIENT , RUN 3 
£=0.0 MEAN SQUARE =0.0727 
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FIGURE 29. PSD OF UNSTEADY PRESSURE COEFRCIENT , RUN 3 
^ = 0.5 LOWER SURFACE MEAN SQUARE =00116 
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£ = 0.125 LOWER SURFACE MEAN SQUARE = 0.0104 



Mean Square x 10 



□ Upper Surface 
A Lower Surfoces 




FIGURE 33. MEAN SQUARE OF UNSTEADY PRESSURE 

COEFFICIENT ON AIRFOIL 

RUN 3 oc = 10° v * 100 ft / sec STEADY FLOW 
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Mean Square x K) 




FIGURE 34. MEAN SQUARE OF UNSTEADY PRESSURE 
COEFFICIENTS ON AIRFOIL 

RUN 4 Oc = 15° \7w 100 ft /sec STEADY FLOW 
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TYPICAL UN FILTERED PSD, RUN 4, £ =0.0 
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Figure 36. Selected Oscilloscope Traces, Run 5 
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TAP 23 




TAP 17 
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Figure 37. Selected Oscilloscope Traces, Run 5 



77 









£ 

i- 

<x> 

X 



>\ 




00 

ro 

id 

u. 



78 



PSD OF UNSTEADY VELOCITY, RUN 5 
1/ = 102.66 V* 38.43 
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FIGURE 39. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 
£=0.0 UPPER SURFACE 




o 

CVJ 



o 

CO 



T 

2 



8 



I r 

8 Q 

aJDnbs UD®w 



5 






80 



Frequency , Hertz 

FIGURE 40. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 
6=0.025 UPPER SURFACE 
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FIGURE 41. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 
^ 0.075 UPPER SURFACE 
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FIGURE 43. PSD OF UNSTEADY PRESSURE COEFFICIENT ,RUN 5 
£.=0.175 UPPER SURFACE 
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Frequency , Hertz 

FIGURE 44. PSD OF UNSTEADY PRESSURE COEFFICIENT RUN 5 
& =0 225 
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FIGURE 46 PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 
6=0.350 UPPER SURFACE 
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FIGURE 47. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 
6=0-400 UPPER SURFACE 
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FIGURE 48. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 
0.450 UPPER SURFACE 
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FIGURE 50. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 
6*0.600 UPPER SURFACE 
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Frequency .Hertz 

FIGURE 51. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 
= 0.700 UPPER SURFACE 
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RGURE 52. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 
t=0.750 UPPER SURFACE 
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Frequency, Hertz 

FIGURE 53. PSD OF UNSTEADY PRESSURE COEFFICIENT , RUN 5 
t=0.0 LOWER SURFACE 
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FIGURE 54. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 
'(j,- 0.050 LOWER SURFACE 
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FIGURE 55. PSD OF UNSTEADY PRESSURE COEFFICENT, RUN 5 
£,=0.100 LOWER SURFACE 
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FIGURE 56. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 
£,=0150 LCWER SURFACE 
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FIGURE 59. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 

6=0. 350 LOWER SURFACE 
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FIGURE 58. PSD OF UNSTEADY PRESSURE COEFFICIENT , RUN 
£•=0.250 LOWER SURFACE 
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FIGURE 59. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 

<5= 0.350 LOWER SURFACE 
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FIGURE 60. PSD OF UNSTEADY PRESSURE COEFFICIENT , RUN ! 
&= 0.450 lower SURF/CE 
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FIGURE 61. PSD OF UNSTEADY PRESSURE COEFFICIENT , RUN 5 
£=0-600 LOWER SURFACE 
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FIGURE 62. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 
& = 0750 LOWE R SURFACE 
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FIGURE 63- PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 
&= 0.850 LOWER SURFACE 
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FIGURE 64. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 5 
**=0.950 LOWER SURFACE 
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FIGURE 67. PSD OF UNSTEADY NORMAL FORCE COEFFICIENT , RUN 5 

MEAN SQUARE =0.100 
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FIGURE 68. MEAN PRESSURE DISTRIBUTION ON AIRFOIL 
<x=!5° fo = 94Hz £,=0085 




tf 0 =0.563 




tf 0 = 0.688 




tf 0 =0.813 




tf 0 =0.938 





tf 0 =l 000 



Figure 69. Time Dependent Pressure, Run 5 
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Figure 72. Pictures, Run 7, Pressure, Velocity 
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Figure 73. Pictures, Run 7, Pressure, Velocity 
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FIGURE 74 . PSD OF UNSTEADY VELOCITY , RUN 7 
U = 100.46 U 2 =25.99 
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FIGURE 75 PSD OF UNSTEADY PRESSURE COEFFICIENT, Rl 
5 = 00 UPPER SURFACE 
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FIGURE 76. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 7 
t = 0.\00 UPPER SURFACE 
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FIGURE 77. PSD OF UNSTEADY PRESSURE COEFFICIENT , RUN 7 
£ =0.500 UPPER SURFACE 
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FIGURE 79 PSD OF UNSTEADY PRESSURE COEFFICIENT ^ RUN 7 
t = 0. 100 LOWER SURFACE 
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Frequency .Hertz 

FIGURE 80. PSD OF UNSTEADY PRESSURE COEFFICIENT, RUN 
£>=0 ' 500 LOWER SURFACE 
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FIGURE 82. TIME DEPENDENT VELOCITY, RUN 7 
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NORMAL FORCE , RUN 7 
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MEAN SQUARE =0.075 
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MEAN SQUARE =0.010 




i 

126 






FIGURE -86. MEAN PRESSURE DISTRIBUTION ON AIRFOIL 
<*? 15° fo = 1 1 Hz t= 0.072 



FIGURE 87. MEAN NORMAL FORCE COEFFICIENT VS. FREQUENCY 
o= 15° u=IOOft/£ec £*0.08 




FIGURE 88. MEAN CENTER OF PRESSURE VS. FREQUENCY 
cx = 15° xf =100 ft/sec £*0-08 
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APPENDIX A 



TRANSDUCER CALIBRATION 

Dyanmic calibration was carried out using the procedure developed 
by Johnson in Reference 14. The results of his work were used to select 
the inside diameter of the tubing and tube length for this investigation. 
The stainless steel tubing used in the model was 0.0625 inches in outside 
diameter, 0.047 inches in inside diameter and 24 inches long. The plastic 
tubing attached to the transducer housing was 0.49 inches in inside dia- 
meter and 11.5 inches long. 

Both transducers were first mounted in a calibration chamber and 
statically calibrated. Figure 89 shows a schematic diagram of the static 
calibration instrumentation . 

A sound driver was then mounted on the calibration chamber and driven 
at known frequencies by the amplified output of an audio oscillator. The 
spectrum from 0 to 500 hertz was scanned to see if any difference in phase 
or gain could be detected between the two transducers. No difference in 
phasing was detected and the ratio of the gains was constant for the 
observed spectrum. 

The next step was to remove the transducer to be used for the investi- 
gation from the chamber, install it in its housing and reconnect it to the 
chamber through steel and plastic tubing identical in length and diameter 
to that to be used in the experiment. The values of phase difference 
and wave analyzer readings for each transducer were then recorded at 5 
hertz increments from 0 to 400 hertz. The gain of the remote transducer 
was then taken as the product of the ratio of static gains and the ratio 
of the wave analyzer readings. Hie amplifiers used to amplify the two 
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transducer outputs were switched and the procedure was repeated. No 
differences in results were noted. Figure 90 shows a schematic of the 
dynamic ca libra t ion ins trumentat ion . 

The spectral interval subsequently used in the analysis of results 
was 1 hertz. A linear interpolation was used to obtain intermediate 
points in the frequency range 0 to 400 hertz and a simple linear extrapo- 
lation made from 400 to 512 hertz. The results are shown in Fig. 10 and 
11 . 
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D.C. VOLTMETER 



FIGURE 89. STATIC CALIBRATION INSTRUMENTATION 
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FIGURE 90. DYNAMIC CALIBRATION INSTRUMENTATION 



APPENDIX B 



ANALOGUE TO DIGITAL CONVERSION 

The available digitizing equipment was not capable of digitizing 
more than one analog signal at a time. This precluded simultaneous 
digitization of the velocity and pressure analogues associated with a 
given tap and ultimately required photographic measurement of phasing. 

The procedure used was based on an available digital computer pro- 
gram. The program takes binary words from one channel of the analogue to 
digital interface of the hybrid and stores them in a given buffer until 
the buffer is full. The contents of that buffer are then written on 
magnetic tape while a second buffer is being filled. When the second 
buffer is full it is written on magnetic tape and the first buffer takes 
on the storage function. The cycle repeats itself for the sample length 
desired. Sample length is expressed in terms of the number of blocks of 
1024 samples per record. The 1024 samples correspond to the storage 
capacity of the buffers. 

The analogue signal to be digitized is taken by a designated trunk 
line fromthe analogue side of the hybrid to the analog to digital con- 
verter. The sampling interval is controlled by interupts into a trunk 
line on the logic board of the analogue. 

For this investigation a sampling frequency of 1024 samples per 
second was chosen in order to simplify the computer analysis of the 
data. For minimum scale error in the conversion process it was desirable 
that the digitized signal be as large as possible. A compromise between 
the desired amplitude of the digitized signal and numerical convenience 
was worked out by making the filter gain be of the order of 100. 
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The filter was used to improve signal to noise-plus -s igna 1 ratios of 
the analogue data* Preliminary investigations had shown that some high 
frequency noise was picked up in the data recording process from both 
the transducer and the tape recorder. ("High frequency' 1 for the purpose 
of this investigation meaning frequencies greater than 500 hertz.) 
Filtering prevented this high frequency noise from distorting the 
resulting spectral estimates. 

For the steady flow runs it was not necessary to preserve any phase 
informat ion and the analogue circuitry used is shown schematically in Fig. 

91. The conversion process in this case was commenced when the "Master 
Clock" on the analogue computer was switched from "Reset" to "Run" with the 
analogue computer in the "Compute" mode. In more specific terms, switching 
the master clock to "Run" enabled the interrupts to the digital computer. 

The unsteady flow runs were digitized using the circuitry shown in 
Fig. 92. A 1000 hertz sine wave was recorded on a third track of the 
data tape. This signal was used to start the interrupts to the digital 
computer at the point on the tape corresponding to the start of the 1000 
hertz signal. The idea was that the velocity and pressure analogues corres- 
ponding to a given pressure tap would be digitized starting at the same 
point on the tape. The large gain on the 1000 hertz signal was used to 
get maximum rise time to the comparator. Utilizing the circuitry shown 
in Fig. 92 also allowed control of the digitizing process from the 
digital side of the hybrid. 

The procedure described above was not adequate to preserve the 
phase information between the velocity and pressure analogue signals, 
Possible reasons include variabilities in the response times of the 
comparator and/or the delay flop in the logic circuitry and/or recording 
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data at too slow a tape speed to be able to get the rise time necessary 
to insure adequate responses of the above analogue devices. 

A method that was not tried would be to use the 1000 hertz signal 
to put the analogue into f, Compute ,f rather than utilizing the delay flop 
and comparator circuitry shown in Fig. 92. This would eliminate at 
least one of the elements in the circuit and possibly both. 

The hybrid has the capability of convert ing 32 analogue channels with 
only a slight time skew. The conversion process is limited by buffer 
size and the speed at which the filled buffers can be written on magnetic 
tape. The time it takes to fill a buffer cannot be less than the time it 
takes to write the contents of a buffer on magnetic tape. This limits the 
sampling frequency. 

In order to sample more than one track at a time the digital computer 
program would have had to be written to store the sampled analogue signals 
alternately in the buffers. The author did not have the time available 
to rewrite the program and the Electrical Engineering Department Computer 
Laboratory does not have the programming support available to rewrite it. 
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CIRCUITRY C ANALOGUE IN COMPUTE] 
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FIGURE 92. LOGIC AND ANALOGUE CIRCUITRY FOR UNSTEADY 

FLOW CONVERSION 





APPENDIX C 



NUMERICAL ANALYSIS OF DATA 



A. ASSUMPTIONS 

The digital analysis was based on the assumption that the pressure 
at any point on the surface of the airfoil was a parametric function of 
the free stream velocity, i.e., 

P = P(x,y,U), U = U (t ) 

The oscilloscope traces shown in Fig. 36 and Fig. 72 indicate the validity 
of this assumption which allowed the establishment of a common time base 
for each pressure analogue sample. 



B. MEAN AND MEAN SQUARE CALCULATIONS 
The mean of a digital array 
re la t ion : 



{ v „} • - - 1 > 



N was calculated from the 



V = 



N 

E V 
n=l r 



where corresponds physically to a quantity measured at time t^ and N 
is the number of samples in the array. For spectral calculations, it was 
desirable that the arrays subjected to spectral analysis have a zero mean. 
This prevents distortion at the low frequency end of the spectrum [Ref. 13]. 
Consequently, the mean was subtracted from the initial array <V > to form 

[3 V * 

a new digital array <U V. This new array corresponded to the unsteady 
part of the analogue signal. 

The mean square of was calculated from the relation: 

N 



u 2 - i i u 2 

N . n 
n=l 
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C. SPECTRAL CALCULATIONS 



1 . Fourier Transforms 

All sampling for this investigation was done at the rate of 1024 
samples per second. Spectral analysis was performed on arrays of 1024 
samples each. This procedure yields a Nyquist Cutoff Frequency of 512 
hertz and an equivalent spectral bandwidth of 1 hertz [Ref. 13]. 

The complex Fourier Transform of the digital array under con- 
sideration was calculated using the relation: 

N-l . nK 

u„ - i E U e^ ni N 
K N ~ n+ 1 
n=0 



where : 

K ’ F C 

= spectral component at frequency ^y^ , K = 1 , ... N/2 

i = V-T 
e = 2.7183. . . 

The inverse transfer functions of the tubing and filter were 
applied to the pressure analogues in complex form in order to get the 
Fourier Transform of the pressure at the surface of the airfoil, i.e., 

p 

K g k“k 

The Fourier Transform of the unsteady velocity was found by 
simply dividing the Fourier Transform of the velocity analogue by the 
transfer function of the filter. 

The inverses of the pressure and velocity transforms were taken 
using the same subroutine as was used to take the transform of the 
original analogue array. Details of the capabilities of this sub- 
routine, called "HARM" in the data analysis program, are given by 
Cooley and Tukey in Reference 15. 
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2 . Power Spectral Densities 




The power spectral density of 
are related by: 




2 V\* 




and the Fourier Transform of 



The power spectral density of a function describes the general frequency 
composition of the function in terms of the spectral density of its mean 
square. Thus the mean square of was calculated from: 

00 



? 2 



= S G. 
K=0 



Af 



K 



where the equivalent spectral bandwidth, Af , was 1 hertz for this 

D 

investigation. Reference 13 has a complete discussion of the properties 
of the power spectral density and its uses. 



D. NORMAL FORCE AND CENTER OF PRESSURE 

The integral relations given in Chapter IV, D, were approximated by 
a simple trapezoidal summation. For example, the mean normal force 
coefficient, defined as 

1 

V/ [C p - c p (§,-h ')] df 

0 



was approximated by: 



C 



N 



23 [AC (-+1) + AC(n)] 

— 2 — 



where : 



AC p (n) = C p (n) - C p (n). 

u L 
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APPENDIX D 



EVALUATION OF EXPERIMENTAL TECHNIQUES 

A. PRESSURE MEASURING SYSTEM 

Dynamic calibration of the transducer and pressure tap leads was 
accomplished using a sound driver to produce reference pressure fluctu- 
ations, (Appendix A). The lower frequency bound on the flat portion of 
the dynamic response curve of the sound driver was 160 hertz with 
detectable response of decreasing amplitude to 25 hertz. 

Below 25 hertz, both the gain of the amplifier used to drive the 
sound driver and the gain of the sound driver were very small. The 
power required to produce any output response from the sound driver 
below 25 hertz caused saturation of the amplifier with the result that 
the output of the sound driver was a distorted signal rather than a 
simple sine wave. The combination of a complex signal, decreasing 
interval between harmonics, and low gain made the wave analyzer more 
difficult to read and the resulting gains subject to greater uncertainty 
than gains measured at higher frequencies. 

Defining error in terms of a maximum percentage of deviation from a 
true value, the estimated error in the system transfer function was 5 
per cent above 25 hertz and 10 per cent below 25 hertz. 

B. TAPE RECORDER AND ANALOGUE TO DIGITAL CONVERSION 

Table II gives examples of the differences between measured recorded 
voltages and their digital equivalents. The differences were considered 
to be insignificant. 
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C. SPECTRAL ANALYSIS 



The power spectral density estimates of the unsteady pressures and 
velocities are considered by the author to be reasonable estimates to 
their true power spectral densities. A detailed discussion of error 
criteria in spectral calculations is given by Bendat and Piersol in 
Reference 13. 

D. NORMAL FORCE AND CENTER OF PRESSURE 

The inability to preserve phase information in the analogue to 
digital conversion process and subsequent measurement of time differ- 
ences from oscilloscope traces made error in the normal force calculations 
difficult to evaluate. 

It was estimated that phasing between the pressure and velocity 
analogues associated with a given tap could be measured to within 1/2 
sample time increment or approximately 17 degrees at 94 hertz and 2 
degrees at 11 hertz. The quality of these measurements would establish 
minimum error bounds on the normal force and center of pressure calculations. 

E. EXPERIMENTAL ENVIRONMENT 

The variation in the ratio of the unsteady to steady components of 
the velocity associated with individual taps for Runs 5 and 7 are given 
in Tables IV and V. The pressure coefficients associated with a given 
tap were calculated using the measured mean velocity associated with that 
tap in an effort to minimize the effect of the variations on the normal 
force and center of pressure calculations. These variations in free 
stream velocity apparently were due to fluctuations in line voltage and 
atmospheric disturbances. 
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TABLE -W. 

RATIO OF SQRT<2)*SQRT(VELOCITY MEAN SQUARE) TO MEAN VELOCITY 



K UIN 

TAP 


LOCATION 


UPPER 


LOWER 


1 


0.0 


0. 9R8 


0 .078 


2 


0.025 


0.991 


0.081 


3 


0.050 


0.090 


0.080 


4 


0.075 


0. *>89 


0.C80 


5 


0. 100 


0.089 


0.0 80 


6 


0.125 


0.089 


0 .0 80 


7 


0.150 


0.086 


C.0«3 


8 


0.175 


0. 090 


C .C83 


9 


0.200 


0.088 


C .080 


ir 


0.225 


0.990 


9.079 


n 


C.250 


0.088 


0.080 


12 


0.300 


C . °87 


0.081 


13 


0.350 


0.091 


f . C 8 1 


14 


0 .400 


0.088 


0.081 


15 


0.450 


0. 989 


r .082 


16 


0.500 


0.991 


0.0 82 


17 


0.600 


0.087 


0 .081 


18 


0. 700 


0.091 


0.084 


1« 


0.750 


0.090 


C .086 


20 


0.8-9 0 


0.087 


0.0 34 


21 


C. 850 


0.986 


0.086 


22 


0.900 


0 .989 


0.0 86 


23 


C. 950 


0.087 


9.0 87 


24 


C.250 


0.088 


0.986 


25 


0.500 


C .090 


0 .086 


THF 


AVERAGE MEAN 


VELOCITY FOR RUN 5 


WAS 


102. 66 






THE 


AVERAGE MEAN 


SQUARE WAS 


38.43 


THE 


RATIO OF UNSTEADY TO STEADY WAS 0 


THE 


AVFRAGE Q WAS 


12.95 
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PATIO OF 
PUN ^ 



TABLE ST 

SGPT (2)*S0RT (VELOCITY MEAN SQUARE) 



TO MEAN VELOCITY 



1 


0. A 


0.069 


C .073 


2 


0.025 


0. 072 


0.0 71 


3 


0.05C 


0.072 


0 .073 


4 


0.0 75 


0.071 


0.072 


5 


0. IOC 


0.070 


C . 0 7 3 


6 


0.125 


0.072 


0 .0 71 


7 


0.150 


0.172 


0 .071 


P 


C. 1 75 


0.074 


0 .0 72 


9 


0.200 


0.071 


0.073 


1C 


0.225 


0.072 


0.071 


11 


0.250 


0.0 71 


0.071 


12 


C. 300 


0.070 


0 .071 


13 


0.350 


0.071 


1.074 


14 


0.400 


0.073 


o .f 71 


15 


0.450 


0.071 


^ .073 


16 


0.500 


0.071 


0.071 


17 


0.600 


0.074 


0 .071 


18 


0.700 


0.072 


0.071 


1 .9 


0.750 


0.069 


0.0 71 


20 


0.810 


0.0 70 


0 .0 7 2 


21 


0.850 


0.173 


O.C 73 


22 


O.o 00 


0.071 


0.0 74 


23 


C. 950 


0.0 70 


0.0 72 


24 


0.250 


0.071 


0.072 


25 


C.50C 


0.0 72 


0 .074 



THE 

WAS 

THE 



AVERAGE 
ICO .46 
AVERAGE 



MEAN VELOCITY FOR RUN 



MEAN 



SGUARE WAS 25.99 
THF PATIO OF UNSTEADY TO STEADY WAS 0 
THE AVERAGE Q WAS 11.63 



0 7 1 R 
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AM EXPER IMFNTAL INVESTIGATION 
INTO THE 

UNSTEADY NORMAL FORCE 
ON AN 

AIRFOIL IN AN OSCILLATING FREE STREAM 
PROGRAMMED BY M. R. BANNING, CAPT. IJSMC. AUG-SEP. 1969 



C 



C 

C 



C 



C 



C 

C 



DIMENSION P ( 4C 9 6 ) ,tl(4C96) , CN ( 1C 24 > , BMP ( 2 , 2 5 ) , 

*BMSP ( 2,25) ,BMU(2, 25) ,BMSU( 2,25) ,DATA( 1024) f S (256) , 

* I NV ( 25 6 ) • M ( 3 ) • Y ( 50),ETA(25),GX( 512) ,X( 25) , KPHA SE ( 2 , 2 5 ) 
*FPFQ(512),TIMEIT(512) ,RANGE(4) 

COMPLEXXB A( 10 24, 1. 1 ) ,G( 512) ,H( 512 ) ,CON 
INTEGER SI DE, HOLE, $F,FRUN, ALPHA 



REAL *8 IT I LI (12)/ 

•* P U N t 

RF AL*8 ITI L2( 12) / 
*PUN 6 

R E At *8 ITI L 3 ( 12)/ 
*«?UN 6 

RFAL*B I T I L4 ( 12)/ 

*'■ k U N 6 

«1AL*8 I T I L5 ( 1 2 ) / 
* RUN 6 

RCAL*8 I T I L6 ( 1 ? ) / 

*RUN 6 

REAL *8 I TILTH 2)/ 
*RUN 6 

PEAL *8 IT I LB (12)/ 
UN (■ 

REAL LABEL / 4H 
REAL*fl LABE1 /OH 
REAL L ABC ( 53) 



•UNSTEADY PRESSURE COEFFICIENTS VS TI 

M.R. BANNING’/ 

•UNSTEACY VELOCITY VS TIME 

M. R. BANNING* / 

•PSD OF UNSTEADY PRESSURE 

M.R. BANNING* / 

•PSD OF UNSTEADY VELOCITY 

M.R. BANNING*/ 

•MEAN PRESSURE COEFFICIENTS VS ETA 

M.R .BANNING’ / 

•PSD OF UNSTEADY NORMAL FCRCF 

M.R. BANNING* / 

•UNSTEADY NORMAL FORCE VS TIMF 

M.R. BANNING* / 

•PSD OF UNSTEADY CENTER OF PRESSURE 

M.R. BANNING’/ 

/ 

/ 



DEFINE FILE 1< 51200, 2, U, IB ) 

READ (5 ,5C 1 ) NR UN, NBSKI P, IFFIL , FR UN .ALPHA 
PE AC (5, 50 3 ) P SLOP E, UMVOLT, GPU, TFMP.P ATM 
REACH 5,502) (FT A( I ) ,1 = 1,25) 

PEAD(5,5C4) (G( I ) , 1=1, 512) 

READ (5 , 5 C 4 ) ( H( I ) , 1 = 1,512) 

READ (5, 505) (LARC< I ), I = 1,50) 

G A SC CN = 2 1 1 6 . 0 / ( 23 769*518. 638 ) 

TLMPR=TDMP+453.688 

RH0=PATM*7C . 7 3 / ( G ASCON*T EM PR ) 

UMEAN=SGRT < 2 .0 *CE 0*5 . 2 32/ ( 2 . 54 *RHO ) ) 



W P IT F ( 6 , 7C 1) NR UN ,NBSKIP,PSLU PE, UMVOLT ,UMEAN 
701 FfRMATdHl, / / / , T 3 2 , * I N PU T PARAMETERS FOR RUN*, 13,/, 

* T 3 2 , * NO. OF BLOCKS SKIPPED BEFORE READING FIRST*, 

*• RF CCRD= * ,14,/, 

*T32» 'TRANSDUCER SL C PE = * , F 5 . 3 , • VOLTS/IN. OF H20*,/, 

* T 3 2 , 'HOT WIRE OUTPUT OF ',F3.1,' VOLTS CORRESPONDS TO' 
* / , T 3 5 , 'A VELOCITY OF • , F6 . 2 , * FEET/SEC') 

WPITF( 6, 7C2) FRUN, TFMP , P A TM , RHC , UME AN 
7 C 2 FORMAT (T32, • FREE STR E AM OSC I L L AT 1 0N= • , 1 4 , • HERTZ*,/, 
*T32, 'ATMOSPHERIC T EMPE RATURE= ’ , F5 . 1 , ’ DEGREES F',/, 

+ T32, 'ATMOSPHER IC PRES SURE = • , F6 . 2 , ' IN. OF HG',/, 

*T32, 'ATMOSPHERIC 0 ENS I T Y= • , F 8 . 5 , ' SLUGS/FT **3 • , / , 

* T 3 2 » 'MFAN VE LUC I T Y = • , F 7 . 2 , • F T/SEC • ) 

WR I T F ( 6 , 70 3) ALPHA 

703 EL RMAT (T32 ,* ANGLE CF ATT ACK= • , 1 3 , • DEGREES') 



REWIND 4 



NCBLK S = ^ 

K I N= 14C 
I r SE T = - 1 
NTOT= 1C?4 
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NDIR=23 
I F F I L = r 
K K = 10 
R( 1 )=KK 
R< 2) = 0 
m ( 3) =>: 

FACl=5.2C2/( . 5 *PHO*PSLOPE ) 

FAC? = U*FAN/IJMVCLT 
FAC3=CABS( H< 1 ) ) 

F AC4*FAC1/ FAC 3 
F AC5 =f AC. 2 / F AC 3 
FAC6=FAC2/ (CAR S( HI FRUN+1 ) ) ) 

C 

PANCE ( 1 ) =0 . 35 
RANG! ( 2)=c .3 
RANGE! 3 ) = 1 5.0 
RANGE ( A ) =- 15. 3 
C 

TC 1C<~2 1=1,512 
FR EQ ( 1 ) = I- 1 

1 002 TIMFITU )=(I-l) = l. C/r 24. 

DC 1 1=1 ,NBSKI P 

1 CALL INUATAIDATA) 

L = C 
C 

DC 11 SIDE =1,2 
DC 11 HOLE =1,2 5 
L = L* 1 
C 

CC 3 I BL K= 1 , 4 
J= ( I PLK-1 ) *1024 
CALL I NDAT A( DATA ) 

DC 3 J J=l, 1024 
1=J+ JJ 

3 U( I ) =0 AT A I JJ ) 

C 

CALL BMC AMU, NCBLKS,RMU(S IDE,HQLF),BMSU(S ICF, HOLE ) ) 
C 

CALL BZERC ( IJ , K IN, K OUT ) 

C 

K PHASE! SIDE, HOLE) =KCUT 
KSTART = KC*UT 
KLAST=KSTART + 1^23 
C 

BMJ( S ICE,HCLE )=FAC5*RMU( S IOE, HOLE) 

B F SU ( S IDE , HOLE ) =( FAC6**2 ) *BMSIJ ( S l I3E, HOLE ) 

C 

CC 0 K=KST ART, KLAST 
I=K-KSTART+1 
6 A( I, 1, 1 ) = F AC, 2*U ( K ) 

C 

If SF T= - 1 

CALL HARR ( A,M, INV , S, ! F SET , IF ERR ) 

C 

A( 513, 1 ,1 ) =A( 513, 1 ,1 )/H(512 ) 

A ( 1, 1 , 1 ) = A ( 1 , 1 , 1 ) /HI 1 ) 

GX(l) = 2. r *A( 1,1,1 ) *CL)N JG ( A ( 1,1,1)) 

PC 21 1=2,512 
KCONJG=NTOT- I f 2 
A ( I , 1 , 1 ) = A ( I , 1 , 1 ) / H ( I) 

CCN=CCNJG( A ( I , 1 , l ) ) 

A ( KC( M JG, 1 , 1 ) = CON 
21 GX(I)=2.0*A(l,l,l) *CON 
C 

SUM=( ,C 
DC 22 1=2,511 
2 2 SLR = SUR+GX( I ) 

SUR=SUR+.S*GX( 1)+.5*GX(512) 

C 

I F SE T= l 

C ^ LI HARM! A, M, IN V , S , I F SE T , I F ERR ) 

DO 7 1=1,1024 
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c 

c 



r : 



c 

c 

c 

c 



c 

c 

c 

c 



7 U ( I ) = A ( I , 1 , l ) 
to P I T E ( 6 , 122) 

WRITE! 6* K 2) „ _ „ 

to P I T E ( 6 , 1(A) NRUN, SIDE, HOLE 
toKITE! 6, 12C ) BVU! SIDE, HOLE) 

CALI UTPLOT ( TI ME I T ,U, 50, RANGE , 1 ,0 ) 

11 CONTINUE 

to P I T E ( 6 , 1 C C ) 

toP I TF ( 6 * 1C 2 ) _ r 

WP I T f ( fe, If A ) NRUN, SIDE, HOLE 
to P I T E ( 6 , 1 C - ) SUM 
DC 61 LINE =1,1 28 
1X1= (L INE-1) *4+1 
f X2-IX 1 4 ** 

I F ( ( L I NE^EG. 44). OR. (LINE. EQ. 89)) toRITE(6,lC7) 

61 WRITE! t, K 6) ( I X , GX ( IX), I X = I X 1 , 1X2) 

WR ITE ( t, 1 IS) 

WPITFC6, 10 2) 

top ITE (6, 1C 4) ( NPUN, SIDE, HOLE) 

to P I T F ( fc, 120 QMU( S IDE , HOLE ) 

Dr P.CCl L I NE= 1 ,89 
I X 1 = ( L INF- 1 ) *4 + 1 
I X 2= I X 1 + 3 

IFIL.EC.44) WRITE(6,1C7) 

11 toPITFI 6,121)1 IX, U( IX) ,IX=IX1,IX2) 

CO lCCl J= 1 , 5 1 2 

" 1 CALL^DRAto'I is'^FPEJ.GX.O^tLABCIL) , I T I L 4 , 6 C .2 ,0.2, 

CALL D8Ato(3)»tl^EIT*U,0*2»LAREl, ITIL2,C .004,0.0, 
,v*,f ,f , 8, 4, 0, LAST ) 

RANGE! 3)=9.1 
RANGE! 4 > = - .=> 0 

L= 0 

CC 1? 5 IDE = 1,2 
DC 1? HOI E = 1 , 2 5 

L= L+ 1 

DC 2 I BL K= 1 , 4 
J= ( I RLK-1 ) *1C24 
CALL I NPAT A ( DAT A) 

DC 2 J J= 1 , 1124 
I = j+ JJ 

? PCI) =D AT A ( JJ) 

CALL E)MEAN ( P, NOBLKS , 9MP ! S IOE, HOLE ) ,3MSP< SIDE, HOLE ) ) 

CALL E ZERO !P , KIN, KCUT ) 

K S T AR T = KOU T + 1 
KLAST=KSTART + 1023 

IK 4 K=KST ART , KLAST 
I=K-KSTAPT+1 
4 A! I, 1, 1 ) = F AC 7*P ( K ) 

CALL T HAPM( A, M, INV, S, I ESET , IFERR ) 

A(E>13,1,1)=A(513,1 ,1)/(G(512)*H(512) ) 

A < 1 , 1 , 1 > =A ( 1 , 1 , 1 ) / < G ( 1 ) *H ( 1 ) ) 

CX ( 1 )= 2 . r * A! 1 , 1 , 1 ) *CONJG( A ! 1 , 1 , 1 ) ) 
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c 

PC 1° 1=2,512 

A(I,l,l)=A(I,l,l)/(G(I)*H(M) 

KC CNJC=NTCT-I +2 
Cl N = CC’NJG( A ( I , 1 , 1 ) ) 

A ( KC CP JG , 1 ,1 ) = CON 

19 GX(I) = 2.('*A(I,1,1) *CON 
C 

SUM = C . ' 

nr 2 c i=2,5ii 

20 SLM=SUM+GX( I ) 

SGM= . 5 *GX ( 1 ) +. 5*GX (512) + SUM 
C 

I F SET = 1 

CALL HAF’M( A,M, INV, S, ITSET, IFEPR ) 

DC 5 1 =1 , 1 C 2 A 

P( I ) = A( I , l , 1 ) 

5 P( I ) = P ( I )+BMP( SIDE, HOLE) 

C 

PI 12 1 = 1,1 >24 
IP=( L-l K1C24+ I 

12 wmPITE( 1* IB) PCI) 

C 

wf II f ( 6, 1 1 9 I 
URITE(6,1C 1) 

WF IT n b , 1 1; 4 ) OR UN , SIDE , HOLE 
U F I T f ( 6 » 1 2 r ) 3 M P C S IDE, HCL F ) 

DC 8r ?r LI NE = 1 ,39 
IX 1= ( L INF- 1 ) *4+1 
I X2= I X 1 + 3 

IF (LINE. EC. 44) UR I TE ( 6 , 1 07 ) 

V r ‘ ° w f I T E ( 6 , 1 1 4 ) ( IX, ?( IX) , IX=IX1, 1X2) 

C 

D( PPL 2 1 = 1, 10 24 
8 00 2 P( l)=P(I)-BMP(SIDF,HJLE) 

UFITE(6,122) 

UP ITE( 6, K 1) 

UP I T E ( 6 » 1C4) NRUN, SIDE, HOLE 
UR I te ( 6 , 1 2C ) B MR ( S I UE , HOLE ) 

CALI. UTPLC'KTIMEIT ,P , 50 , RANGE , 1,0) 

C 

Uf I TE ( 6, IOC ) 

UP I T F ( 6, 1C 1 ) 

UPITE<6» 1C 4) NRUN, SIDE, HOLE 
UF ITE ( 6, K 5) SUP 
Df f>( L I N E = 1 , 1 2 fi 
IX1=(LINE-1)*4+1 
I X 2= I X 1+ 3 

IF ( (LINE .FC.44 ) .OK. (L INE.EQ.89) ) WR[ TE ( 6, 1CM 
60 UFITE(b,lC 6) ( I X , G X ( IX) , I X = I X 1 , 1X2) 

C 

I'D P.‘ 5 0 1=1,512 
3 '.50 GX ( I ) = GX ( I ) / SUM 

CALL D R A U ( 4 8 C » F R E 0 » G X , C , r , L A B C ( L ) , ITIL3,6C.0,0.2, 
* > ,C ,0 ,C , 8 , 4 , w> , LAST ) 

C 

13 CONTINUE 
C 

UK I T E ( 2 , 50 CO ) RWP 
C 

DC 5C I T = 1 ,1024 
DO 51 HOLE=l,5C 
I r = ( K LF-l ) * 1* 2 4 + I T 
PFADIl'IBI Y(HCLF) 

S I DF = 1 
I F OL F= HOLE 

If (HOLE. GT. 25) l HC LE =HCL E -2 5 
IF (HOI E .GT .25) SIDE = 2 
51 Y ( HOLE) =Y( HOLE )-BMP( SIDE, I HOL E ) 

C 

UP ITF ( 2, 50( 3 ) Y 
C 
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c 

c 



DC 5 2 HOL E - 1 ,25 
Y ( HOLE ) = Y( HUL E + 2 5 ) -Y( HOLE) 

52 X (HOLE )=Y( HOLE)*ETA(HOLE) 

CALL QTFG( ETA, Y,Y,NDI Ml 
CN(IT)=Y(NCIM) 

CALL iJTFG(ETA,X,X,NDIM) 

C A T A ( IT) = y < NO I M) 

50 CENT INUt 

SUM 1 =C .C 
SL M2 =C .0 
DU 27 S I DE =1 .2 

sf Ml=SUMl+eMU( S IDE, HOLE) 

SUM2=SUM2+RMSU(SIDE»HCIF) 

2 7 CONTINUE 

OC 36 S I DE = I ,2 

36 PMSU ts’lDE^ HOLE ) = SORT ( 2.0 ) *SQRT < BMSUt S IDE , HOLE ) )/ 

* BN U SI HE, HOLE ) 

EM UA = SUM 1/50.3 
PM SO A= SUM2 / 50 .0 
Uf A=RMUA 

RAT 10= SORT ( 2 . 0 *UMS A ) / UMA 
QM EAN= ,5*RH0*UWA**2 

pp t o hPLE= L t 2 3 

Y ( HOLE ) = BMP( 2~,HOLE) -BMP! 1 , HOLE ) 

29 X ( IIOL F )=Y( HOLE )*FTA( HOLE ) 

CALL QTFGI ET A ,Y ,v , NCI M ) 

F NME A K=Y ( NDI M) 

CALI QTFGC ETA, X,X,N0I M ) 

CMMF AN=X t NC I M ) 

il C A TM I )=(CMMEAN+OATA< I ) ) / ( FNME AN + C N ( I ) ) 

M P R I K ^ = 1 

CALL BMEAN (DAT A, N06LK S , CPM, CPMS ) 

V> P IT F ( 6 1 Ilf) NRUN, FRUN, RATIO, UMA, ALPHA 
w KITE! 6 111 2) F NME AN 
DO 32 L I N F = 1 » 8 9 
I CM = ( L If E-l ) *A + 1 

IF (l1nE N 1f-0. AA) V.RITEI6, 107) 

32 Wl KITE (6,113)1 I CN,CN(ICN)» I CN- 1 CN1 » ICN2 ) 

W P IT E ( 6 * 1 1 C ) NRUN, FRUN, RATIO, UMA, ALPHA 
UP ITE ( 6, 1 1 1) C PM 
DP 33 L INF =1,39 
IC Fl=< LINF-1 )*A + 1 
IC P2= I CP 1 + 3 

v i[iW«6.ilM(ICP,0«TMlcj!?ICP=ICPl.ICP?) 

WRITE (6, 21 8) NRUN 

37 WPITE(6?2 F 1 A) f HOLE , ETA < HOLE ), TRMSUT S TOE , HOL E) .SIDE -1*2) 
y, P I T E ( 6 * 2 r ' P ) NRUN * U M A , (IMS A * RAT I 0 * QM EAN 
WRITE! 6*213) NRUN 

28 WF I T t ( 6* 2 1 a} ’ hCLE , ET A ( HOL E ) * ( BMP ( S IDE » HOLE ) , 

* S I TE = 1 * 2 ) 

WPlfE(6*215) NRUN, FNM EAN 

DC 3A 1 = 1,1 D2A 
3 A At 1,1,1 ) =C N ( I ) 

CALL T H ARM! A,M, INV, S, I T SET , I FERR > 
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Sl'M=' .<■ 

DC 35 1=1,512 

GX(I)=2."*A(I,1,1) *CONJG( A (1,1,1)) 

35 SI M=SUM*-Gy ( I ) 

Sl> = SUN-. r *GX< 1 ) - . 5 *GX ( 5 1 2 ) 

wFITf(h,l' C) 

WRITE (6, 115) 

WRITE! 6,110 NKUN ,FKUN, RATIO, DMA .ALPHA 
WE I T E ( 6 , 1C 5) SUM 
L'C ?P L I ME =1 ,128 
I C M = ( L I ^ t- - 1 )*4+l 
I C N2 = I CN1 + 3 

I f ( (LINE.FC.44) .OR. (LINE.FQ.3A) ) WRITF(6,107) 

3b WF ITF( C, 1< 6) ( ICN,GX( ICN) , ICN= ICN1 , ICN2 ) 

PC ICC A HOLE = 1 ,23 
1004 Y ( HOLE )=-RMP ( 1 , HOLE ) 

CALI OPAW(23,ETA,Y,l,3,LABEl,ITIL5,3.2,l.C, 
*G,C,C,C,8, 4,0, LAST) 

C( ICC 5 HE LF = 1 ,2 3 
n‘5 Y(HOLE)=-BMP(2 .HOLE) 

CALL DRAW! 23 ,ETA , Y , 3, 5.LABE1 , I TIL5 , O 2 ,1 . 3 , 

,C,D, r ,8,4,0, LAST ) 

CALL DRAW ( 3^ , T IMF I T , C N , O 2 , L AB E 1 , I T I L 7 ,C . 0C4 , 0 . •> , 
, r , C, 0,8,*, 0, LAST ) 

DC ICC 6 1=1,512 
KC6 G V ( I ) = G X ( I ) /SUM 

CALL OP AW ( ABO , FRFQ.GX ,o,C , LABEL , IT IL6 , 6-' .0 , . . 2, 

^ ,C ,C,‘. ,8,5, 3 , L A S T j 

DC 3 C I = 1 , 1 3 2 A 
35 A ( I , 1 , 1 ) =(> AT A ( I ) 

I FSFT=-2 

CALL HAPM(A,M,INV,S,IFSET,I C ERR) 

SUM=r.o 

DC AC. 1=1,512 

GX ( I ) = 2.C<"A( I , 1 ,1 ) *CONJG( A( I , 1 , 1 ) ) 

AO SLM= SUM + GX ( I ) 

SL M= SU M- . 5 *GX ( 1 )-.5*GX(512) 

WF IT F ( 6, 1C<* ) 

WRITF(6,116) 

wRITE(6,110) NRUN, F RUN, RATIO, UM A, ALPHA 
WFITE(6,1C5) SUM 
DC A 1 LI NE = 1 , 1 28 
ICP1=< LINE-1 1*4+1 
IC P2= ICP1 +3 

I F ( ( l INt.FQ.AA).0R.(LINF.EQ.8 Q ) ) WRITE(6,107) 

A 1 WPITF(6,1C6)( ICP,GX( ICPI , I CP = I CP l , ICP2) 

DC ICC 3 1 = 1,512 
no-8 G X ( I ) = G X ( I J/SIJM 

CALL CPAW( A8 j, FRED, GX,0, '.LABEL, I TIL 8,60.0,3. 2, 

*C ,C,C, 0,8, 5,0, LAST) 

DC A 2 1 = 1, I) 2 A 

A 2 P ( I ) =CM I ) 

(< K P = P* 

CALL RhAFM ( P ,KKR , I NV, S , IFERR ) 

WPITF(6,lCO 

WRITE! 6, 1 1 7 ) ( P ( I ) ,1 = 1,1326) 

PC A? 1=1, 1 0 2 4 
A3 P ( I ) =0 AT A ( I ) 

CALI PHARM ( P ,KKR, I NV, S, IFPRR ) 

W f< I T F ( 6 , 1 C C ) 

WRITE(6,117) ( P ( I ) , 1 = 1,1026) 

CALL PhARMOJ.KKR, INV, S, IFERR ) 

W P I T F I 6 , 1 C C ) 

WRITE! 6, 117) (U( I ) , 1 = 1,10 2 6) 
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ro FCPM AT ( 1 ) 1 , ///////. T 53 , ' T A 6 L t , / , 

*T 46 , 'POWER SPECTRAL DENSITY' ) _ icmti , 

pi FORMAT (T 42 »' UNSTEADY PRESSURE COEFFICIENT ) 
FORMAT (TAQ ,' UNSTEADY VELOCITY') 

F 0 PMATU 42 ,'RUN' ,1 2 .T 53 ,' SIDE' , I 2 .T 66 , TAP , 13 ) 
FORMAT (TAB ,' MEAN S QU AR E= ' , F7 . 3 ) 

FORMAT ( lh 1 | / / //// / * T 4 7 , ' TABLE^ * * ( CONT INUED) ' »/) 

Rf ic. 5 NTSM 



I ' 2 
1 0 4 
i :5 
1C 6 
1C 7 
13 8 



110 



FORMAT ( T 12 ,'RUN' , I 2 , 2X , ' FREQUENCY* ' » I 2 » ' HZ^ ,2X , 
..MOAN VELOCITY=',F6.2, ' FT/ScC', 



111 
1 1 2 

I 1 3 

I I A 

115 

116 
117 
115 
120 
121 
122 
2 A 8 



* ' E PS = * ,F5. 3.2X 
* 2X , 'ALPHA*' ,12,' 
FC RM AT (T A2 , ' MEAN 
FC RM AT ( T 3 G t ' MEAN 
FORMAT ( 11 X ,A( • CN ( 
FC RM AT ( 1 1 X , A ( 'CPI 



DEG' 
CENTER 
NORMAL 



t 2 X , ' SF = n2A 



. . . HZ' ) 

OF PRESSURE*' ,F6.3 ) 

FORCE COEFFICIENT* ',F6. 3) 



• ,13,' )=' .F6.3.11X) ) 
|.I8.‘)='»F5.3,11X)) 



F C RM AT ( T A 1 , ' UNSTEADY 



, 13, 

" NCRMAL 
CENTER 



FORMAT (TA3 , ’ UNSTE ACY 
FORMAT (IX, 1 > E 1 C • 2 ) 

FORMAT (1 HI ,///////, T53, 'TABLE' ) 

FC RM AT ( T 5 r , ' MEAN* ' ,F8 . 3 ) . Y . . 

FCRMAT(11X,4('U('»13,')='»F7.3,1QX)) 

FCRM.AT ( 1H 1 ) 



FORCE COEFFICIENT') 
OF PRESSURE') 



THE AVERAGE MEAN VELOCITY FOP RUN', 13,/, 



FORMAT (1 

*T^2»'THE AVERAGE MEAN SQUARE ViAS * , F7 . 2 , / , 

XT 32, 'THE RATIO OF UNSTEADY TO SThADY WAS'»F7.4,/, 
**T?2,'THE AVERAGE Q v>AS'»F7.2). 

21 3 FORMAT! 1 H 1 ,/////// » T 5 3 » ' T ABL " ' ,/ , 

-T32,'MFAN PRESSURE COEFFICIENTS',/, 



* T 3 2 , ' P UN ' , I 3 , / , 



•CP UP°ER ' 



*-132, 'TAP ' , T39, 'LOCATION' , T 5' 

* ' T 6 C , 'CP LOWER') , _ v c , 

214 E[5^T(^ 2 !ffil:‘fftE 5 AliN^RAiC 5 F6Ec6 3 COEFFlCieNT F OP 



215 



RUN* ,13,/ 



~ ~> . 



y *T*2»'IS',F6»3) 

? 1 8 FC PM AT (IH1, ///////, T53, 'TABLE',/, 

AT 32 , 'RATIO OF SORT (2)*S CRT (VELOCITY 
MEAN SQUARE) TO MEAN VELOCI T Y',/, 
*T3 2 » 'RUN' ,13,/, 



*T32* 'TAP' !t39,'LOCATICN' ,15'’ , 'UPPER' ,T6(,' LOWER' ) 



<=,r i FORM AT ( 5 I 5 ) 

5 j 3 FC RM AT (5F1C.O 
50 2 FORMAT ( 1 6F 5. C ) 
52 4 FORMAT (8EIC .3) 
52 5 FORMAT ( ?C A4 ) 
5)30 FORMAT ( 5CA A) 

ST CP 
END 



SLePfUTINE OTFG(X,Y,Z,NDIM) 

DIMENSION X! 1 ) ,Y ( 1 ) ,Z ( 1 ) 

S U M 2 = C . 

DL 2 I = 2 , N D I M 

SUM2*SUM?+0. 5 * ( X( I ) - X ( I- 1 ) ) * ( V ( I ) ♦ Y ( I - 1 ) ) 
Z ( 1-1 ) =SUM 1 
7 (NDIM ) = SUM2 
R F TUP N 
END 
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SIT ROUT I NE HME AN( P, \OBLKS , BM, BMS > 

C/ 1 MENS ICN P ( 1 » 

MK =C 
3MS=C 
BN =/ “ . C 

1 DC b I=l,NGiLKS 
J = ( I -1 1*1 . 24 
SUM=C . C 

nr ? i i = i f i '24 
JJ= J+ I [ 

IF ( M K .FC. II GG TU 2 
SCM= SIJM + P ( J J )**2 
GO TC 3 

2 SUM= SUM + P ( J J ) 

3 CCMIMJF 
SCN=SCM/ir 24 

IF (VK .EC. ') GO TO 4 

BF S = 5MS+ SUN 
GC TC 5 

4 BM = B M + SUN 

5 CONTINUE 

IF (MK .EC. ?) GC TT h 

BNS=PMS/NCBLKS 

G( TC 3 

6 3M = PK/N(.FL KS 
CO 7 I=1,JJ 

7 P( I ) =P( I )-BM 
MK =1 

GC TC 1 
3 CONTINUE 

I"I Ff RM AT ( ?^X, •N')RLKS = • , I 2, /, 27 X, • MEAN=* , FI 7.4,/ ,2C X , 

* = • * F I f .4,//) 

133 FC RMAT (2' X THE FOLLOWING ARRAY LISTS THE INPUT T I MF H 
*,/, 2C X » 'ZERO MEAN AND NORMALIZED WITH RESPECT TG TH 
*E 1 ,//) 

1C 2 FORMAT! IX, ] ">E I 1.3) 

RETURN 

END 



S L PR OU TINT INDATA(CATA) 
OIMFNSICN DAT A ( l r ?4 ) 

1 FORMAT (16! 64A4 ) > 
RCeDU.l) DATA 
return 
END 



St PF CUT INF 3 ZERO ( X, KIN, KOUT ) 

I I MF NS I UN X ( 1 ) 

X 7 EE 0= 1 . r 
KOUT = 1 

rn l k= 2 1 , k i n 

IF (X (K ) .GT . 3. " > GO TO 1 
I F ( X ( K + 1 ) .LT. r.O) GC TO I 

I F ( ( X! K-l ) .GT. O. 1 ' ) .OP. ( X ( K-2 > . GT.C . T ) ) GO TO 1 

CK = APS ( X (K ) > 

CKP1=APS(X (K + l n 
IF (CK.LT.CKPI > GO TO 2 
K T RY = K + 1 
X TR Y = C KP 1 

IF(XTkY.lT.XZERO) GO TO 3 
GC TC 1 

2 KTP Y = K 

X T R Y = C K 

IF (XTRY.LT . X Z EPC ) GO TO 3 
GC. TC I 

3 KOUT =K TR Y 
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F F TRY LTHtL'K X,Y,NDATA, RANGE ,KKZ,MODCUR) 

A4A If (MCDCUR. EQ.O.OR.MOOCUR.EQ.l ) JSET=9 
JSET=JSET+1 
IF(JSET.GT.A) J SE T = 1 
01 7 Of I =1 ,NOATA, KKZ 

I PTX = 6'-.* ( YMAX-Y (I) I/YRANGE + l .5 
I PTV=F. X( I ) - XM I N ) / XRANGE+ 1 .5 
IFdPTX.GT .6 1. OP. IPTY.GT.B1) GO TO 70 
IT ( I PTX. Lt . 9 .OP. I PTY. LE.C I GO TO 7, 

GP I D ( I PT X , I °T Y ) = XCHAR(JSET) 

GL TC 7u" 

7 C IERR=IERP+1 
7;0 Cf NT I MJE 

C 

C COMPUTE PROPER SCALE NUMBERS 

C 

I F ( MCDCUR . E «• 3 • CR . MCDCUR .60.1) GO TO 800 0 
I F (MCDCUR. FQ.2 ) RETURN 
Gf TO 022 

s:cn XINCF-XPANGE/^. 

Y I N( R= YRANGE /6 • 

XSCAl L ( 1 ) = X M A X 

Y SCALE ( 1 ) =YMAX 
DU PC 1-2, 5 

8 Z XSCALFf I ) = XSCALF( I-D-XINCR 

DC PI 1=2,7 

81 YSCALF ( I ) = YSCALF( I-D-YINCR 
C 

C OUTPUT SECTION WITH GRA^H 

C 

I F (MrDCUP. EQ.O .CR .MCDCUR.CQ. 3 )G0 TC 022 
RETURN 

17 F( RM AT ( I2X , 

I IP El ".3,A( K X,E1 A . 3)/15X,2H**,8(lfri+****+**** 

122 W F I T E ( 6 , 1 7 ) XSCALE( 5) , XSCALE ( A) , XSCAIE( 3 I , XSCALEI 2 ) , XS 
11=1 
I =•" 

DC K1 I K = 1,61 
IT (1)^1 ,«1 ,92 

91 WRITE! 6,18) YSCA L E ( I I ) , ( GR 1 0 ( I K , I X ) , I X = 1 , 8 1 ) , YSC AL E (I I 
13 Ff-RMAT (3X, IP 

1 F10.3,2X,IH+,1X,81A1,1X,1H+,2X,F19.3) 

11=11+1 

GC TO K? 

92 W F I T E ( 6 , 1 9 ) < GR ID ( IK , I X ) , I X = 1 , 81 ) 

19 FORMAT ( 15X , 1H* , IX, 81A1 , IX, IH*I 

102 1=1+1 

IF ( 1-12 I 1C 1,10 3,1 93 
193 I=C 
131 CONTINUE 

W F I T E ( 6 , 2 2 ) X SCALE (5) , X SC ALE ( A ) , X SC AL E t 3 I , XSC ALE ( 2 ) , XS 
22 FC RN AT ( 1 EX ,2H**, 8( 10H+*******-** I , ^H + **/ IP 

1 12X ,E K. 3, A( 10X 

IF ( I ERR) 1C0C, 1C 0 : , 1 *' 1 
19' 1 WP IT F( 6,2c > I ERR 

2 rCPNATUCy ' NUMBER OF °C I NTS OUT OF RANGE =' I A ) 

1 : ' 0 R F TURN 
END 
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i r • v ( 1 )=<: 

CL 98r l. = 1 , M T 
I N V ( L V 1E XP+l ) = ^TLEXP 

DC °7r J = 2 fLMIFXP 
J j = J + L P 1 EX F 

0 7 0 I N V ( JJ ) = I N V ( J) + PTLEXP 
PI LEXF =MTLEXP/2 
-NBC LM1EXP=LM1 EX P* 2 
7 3 2 If ( I F bFT) 1 2, 895, 12 
EE C 



SLBROUTINE UTPLCT (X ,Y , NO AT A , k ANGF , KKZ , MQGCUR ) 
Cl MF NS ION GR ID ( M , 81 ) , XSCALFI 5 ) , YSC ALE ( 7) 

DIMENSION X (1),Y (1)»RANGF(4) 

INTFGEK*2 07 ID , BLANK.OCT, XCHAR (4) /1H. ,1H+, 1H*,1HX/ 
DATA TUT , BLANK //.A BA. , / 

G k ID IS THE MATRIX USED TO PLOT THF POINTS 
ICRP=~ 

X»> AX=R ANGE ( I ) 

XM IN=P ANGL ( 2 ) 

YM AX= R ANGF ( j ) 

YN IN=h ANGF ( A ) 

CHECKING X AND V r>CINTS AND ^LOTTING THOS^ OUT CT 
AT THF MARGIN 

DO 3.' 1 = 1 , NO AT A,KKZ 

I M X (I)-XMAX) ?. r 5,225,22 C 
1 2 0 X ( I ) =XMAX 
IE «R = I F RF+ I 
GC TO 21' 

2 :■ 5 I F ( X ( I ) - X M I N ) 2 ~ 3 , 2 I ~ , 2 17 
2*3 X ( I > - X M I N 
I F RR= IEFR+1 

3 1C I F ( Y U)-YMAX)215,215»?12 
312 Y (I ) = YMA X 
I ERR= It’PR + 1 
GOTO ? : 

215 I r c Y ( I J-YMIN >217, 3C,3 r 

217 Y ( I ) = Y M I N 
IE RF = IE KR+ 1 

30 CONTINUE 

PLOTTING X \ND Y AXIS , IF NECESSARY 

XF ANGf = X F A X-XM I N 
YE' ANCC=YMAX-YM IN 

BLANKING OUT MATRIX-(GRIO) 

cc ?"•: i = i,6i 

no ? ci j j = i , 8 i 

0 1 GF I0( I , J J ) ='JLANK 
3 '0 CC NT I NUE 

ytfst=ymax*ymin 

XT EST= XMAX *XM IN 
IF (XTEST) 1 ,222,2 22 

2 2 2 IF ( YTE ST ) 3 33 , AAA , AAA 

I IYAXIS = «E.M-XMIN)/XRANGE+1.5 
1)1 AC 1 = 1,61 
4 0 GP I n ( I , I YA X I S ) =00 T 
GC TCj 222 

3 3 3 I XAX I S=6. . TY^4 X/YRANGE+1 . 5 

Of.' 6f I = ! , P 1 
60 GP 10 ( I XAX I S * I ) =00T 

PLACING POINTS IN THEIR PRO°FP GRIG PUSlTlGNs 



R ANo 
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GC TG r '^C 

920 IP 1= INV( J P 1 ) /NTVN1 
9 30 I=2*( IPPl+IPlJ + l 

IP ( J- I ) 34^ , 9^5 » 343 
340 T = A < I ) 

A ( I ) = A ( J ) 

A ( J ) = T 
T=A( 1 + 1) 

A ( I + 1 ) = A ( J + 1 ) 

A ( J+ 1 ) =T 
345 CGNTINUr 
950 J = J + 2 
36G JJ1=JJ1+JJC1 
C END CF JPP1 AND JP2 

C 

87 0 J J 2 = J J2 + J J C? 

C END GF JPP 2 AND J P 3 LCCPS 

C 

980 J J 3 = JJ3 + JJ03 
C FNC CF JPP2 LOCP 

C 

390 IF ( I PS 6T )891 r 895, 395 

391 l)( 89? I = 1 , NX 

992 A< 2*1 ) = ~ A ( 2 * I > 

395 RFTIJPN 

THE POLLUTING PROGRAM COMPUTES THE SIN AND INV TABLES. 

iro MT=MAX0(M( 1) » M ( 2 ) » N ( 3 ) ) -2 
NT = P AXf ( 2, MT ) 

9C4 IF ( WT-1G) 906 ,906,905 
9C5 IFEFP = 1 
GG TG 695 
906 I FERRIC 
N T = 2 **MT 
NT V2=NT/2 

SPT UP SIN TABLE 
THFTA=PIE/2*ML + 1 > FOR L= 1 
9 1C THFTA=. 7853931 634 

JSTFP=2**< PT-L + 1 ) FOR L = 1 
JS TP P = NT 

JD I F =2 ** ( PT-L) FOR L = 1 
JDIF =?**(MT-L ) FOR L=1 
J [' I F =■ N T V 2 

S< JDI F ) = SIM THETA) 

DO 9 5F L = ? , M T 
THETA=THFT A/2 . 

J S TE P2 = J S T F P 
JSTPP=JDIF 
JG IF = J STEP/2 
S( JDIF) = SIM THETA) 

J C 1=NT-JDI F 
S(JCl)=CrS (THETA) 

JLAST=NT-JSTEP2 
IF (JLAST - JSTEP) 950,920,920 
920 DG 94C J=JSTEP, JLAST , JSTEP 
JC=NT- J 
JI =J +JDI F 

940 S(J0)=S(J)*S(JC1)+S(JDIF)*S(JC) 

950 CONTINUE 

SIT UP INV(J) TABLE 

960 MTLFXP=NTV2 
C 

C PTLFXP=2 ,VA (MT-L) . FCR L = 1 

LM IT XP=1 
C 

C L P 1 F X P = 2** ( L- 1 ) . FCR L = 1 
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C PIT — p F V F F S F 0 . THE FOLLOWING ROUTINE PUTS THEM IN OR OF 

NT SO=NT*NT 
M2MT = M 7 -MT 

iSC ir(M3PT) 3 7j,36F,3fcO 
C 

C 1*3 GP . CP Ej. MT 

36 3 I C C 3 = 1 

N ' J VN T= N3 / N T 
M I NN3=NT 
Gl 'TO 3R3 
C 

C M3 LFSS THAN MT 

370 I C C? =2 
N ? VNT = 1 
N T VN 3= NT /N 3 
PI N N 3 = N 3 

3 30 J J 03 = NT SC/ N 3 

M2MT=M2-MT 

•+ 50 If ( M2MT )47 ; , 46' ,46* 

C 

C M2 GP. GP F]. MT 

460 I G C2 = 1 

N2VNT=N2/NT 
M I NN 2= NT 
GL TC 480 
C 

C M2 LFSS THAN MT 

4 7 C IC.C2 = 2 

N2 VF T= L 
NT VT. 2= NT /N 2 
MI NM2 = N2 

483 J JT 2=N TSU/N? 

Ml MT=M1-MT 

65C IF (*M mt) 570,560, 56C 
C 

C Ml GP. CJK FO. MT 

560 I G C 1 = 1 

N 1 VN T= N 1 / f- T 
M I NN 1 = NT 
GO TO 560 
C 

C Ml LTSS THAN MT 

6 7 C I GC 1 =2 

N 1 VN T - 1 
NT VM = NT/N1 
w I N N 1 = N1 

IPO J J D 1 = NT SC / M 
40 0 J J 3 = 1 
J = 1 

fin ° 8! J PP 3= 1 , N3 VN T 
IPP3=INV( JJ3) 

PC B7C JP?=1,MINN3 
G( ’ TO ( 6 1'. , 62' 1 ) , IGC3 
MC IP3=INV(JP3)*N3VNT 
G( TO 6 3C 

3 2G IP 3= INVC JP 3 > /NTVN3 
630 I 3 = < IPP^+I P3) *N2 

7"o j j 2= l 

DC: 87 ' J PP 2= 1 , N2VNT 
I PP2 = I NVC J J2 > + 13 
()( 6 60 JP 2 = 1, MINN 2 

GC TP (7 K , 72: ) , IGC2 
71C IP2=1NV( JP2I*N2VNT 
GO TO 730 

020 IP?=IMV( JP2I/NTVN2 

7 30 I 2 = ( IPP2+I P2 )*N1 
8C A JJ1=1 

DC ft<- J PP 1 = l , MVNT 
I P P 1 = I N V ( J J 1 )+ 12 
DC 8 5C JP 1 = 1 , M INN 1 
GC TO (Pit ,320 ) , I GC 1 
310 IP 1= INVC JP 1 ) *N1VNT 
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17 0 I ? C = - I 3C 

W3< 1 )=-S( I 3C ) 

W?(2) = S( I3CC> 

gc m 2^r 

180 W3 ( 1 ) =- 1 . 

W 3 ( 2 ) = r . 

GC TC 2 f C 

3* I IN THIRD QUADRANT 
1D0 I 3CCC=NT + I 3CC 
I 3CC = -I3CC 
W3 (1 ) = -S ( I3CCC ) 

W3 ( ? ) =-S ( I 3CC ) 

200 ILAST= IL+JJ 

nr 2 ?r i=jj,ilast,idif 

Kl A S T = K L + I 

C( 22f K = I ,<LAST »2 

K 1 =K +K B I T 

K 2 =K 1+ KP I T 

K3=K 2 + KBIT 

PC TWO STEPS WITH J NOT 0 
A(K)=A(K)+A(K2)*W2 
A(K2)=A(K)-A(K2)**2 
A ( K 1 ) = A ( K 1 ) * W + A ( K 3 ) * W 3 
A ( K3 )= A ( Kl )*W- A( K3 ) *W? 

A ( K) =A (K ) + A( Kl ) 

A C K 1 )= A ( K ) — A ( K 1 ) 

A(K2)=A(K2)+A(K3)*I 
A ( K3 ) = A( K2 ) — A ( K3 ) *1 

P=A(K2)*W2 11 ) - A ( K 2 + 1 ) * W 2 ( 2 ) 

T=A(K2)*W2(2) + A(K2+1)*W2( 1) 

A(K2 ) = A ( K ) -R 
A ( K ) =A ( K ) + R 
AIK2+1 ) = A ( K+ 1 ) - T 
A ( K+ 1 ) = A ( K + 1 ) +T 

R=A(K2)*W3(1)-A(K3+1)*W3(2) 

T=A(K3) y 'W3 (2)+A(K3+l)*W3( 1 ) 

A W P = A ( K 1 ) * W ( 1 ) — A ( Kl + 1 ) * W ( 2 ) 

AWI=A(K1KW(2)+A(K1 + 1)*W(1) 

A ( K? ) = AW R-R 
A( K 3+ 1 ) = A W I - T 
AIK! ) = AWP + R 
A ( Kl + 1 ) = A W I + T 
T= A ( K 1 ) 

A( Kl ) = A(K ) -T 
A (K) =A (K )+T 
T=A( Kl+1) 

A( Kl+1 ) = A( K + l ) - T 
A ( K+l ) = A ( K + l ) + T 
R = -A ( K 3+ 1 ) 

T = A ( K3 ) 

A( K3 )=A( K2 )-R 
A( K2 1 = A ( K2 )+R 
A( K3+1 )=A( K2+1 )-T 
720 A ( K 2 + 1 ) = A ( K2+1 )+T 

Eh C OF I AND K LOOPS 

23G J J = J jn i r + j j 

FNC OF J-L00P 

235 JL AST=A+JLAST+3 
2A0 C( NT 1 MIF 

END OF L LOOP 

250 CC NT INUC 

Eh D OF IP LOOP 

WE NOW HAVE THE COMPLEX FOURIER SUMS BUT THEIR ADDRESS 
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R - — A ( K 2 + 1 ) 

T = A ( K 2 ) 

A ( K 2 ) = A( K ) -P 
AIK) = A ( K ) f R 
A ( K2 + 1 )=A( K+ 1 ) -T 
AIK*! )=A(K4l ) + T 
C 

AWR=A( K1 )-A( Kl+1 ) 

AW I = AIK1 + 1 )+AIKl ) 

R = -A(K3)-AIK3 + 1 ) 

T = A ( K 3 ) - A I K 3 1 ) 

A (K? )= ( AWR-R ) /R00T2 
A ( K 3+ 1 ) - ( A W I - T ) / R C C T 2 
A ( K 1 ) = < AWR+R ) /ROOT 2 
AtKl+1 )=( AWI+T ) / R 0 0 T 2 
T = A ( K 1 ) 

A ( K 1 ) = A ( K )-T 
A ( K ) = A ( K ) + T 
T = A ( K1 41 ) 

A ( K1 4 1 ) = A ( K4 1 )-T 
At K + 1 ) = A ( K 4 1 )4T 
R = -A ( K 3+ 1 ) 

T= A( Y ? ) 

A ( K 3 ) = A ( K? )-P 
A(K2)=A(K2)4R 
A ( K34l ) = A ( K ?4 1 ) - T 
*5 A ( K241 ) = A ( K 2+ 1 )4T 

IF(JLAST-l) 235, 235,^2 
90 JJ= JJ 4 JJ9IF 
C 

C NfW DC THF REMAINING J'S 

DC 23C J = 2 »JLAST 

FETCH W • S 

OFF- W=W**INV(J), W2=W**2 , W3=W**3 
96 I = I N V ( J 4 1 ) 

93 IC=NT— I 
wm=S( i c; ) 

W ( 2 ) =S ( I ) 

12=2*1 
1 2 C = NT - I 2 
IE ( I 2C ) 12E , I 19 , 1 
C 

C 2* I IS IN FIRST QUADRANT 

l W2 m = S (1 20 
W 2 ( 2 ) = S ( 12 ) 

GC TC 13E 
1 IC W 2 ( 1 ) = , . 

W 2 ( 2 ) = 1 . 

GC TP 1 3 C 
C 

C ?* I IS IN SFCOND QUADRANT 

120 I 2 CC = I 2 C ♦ N T 
I 2C =- 1 2C 
W 2 ( 1 ) = -S ( I 2C ) 

W 2 I 2 ) = S ( 1 2 CC ) 

139 13=1+1 2 

I 3 C = NT - I 3 

IF I I 3C ) 16C , 150 ,1A ? 

C 

C 13 IN FIFST QUADRANT 

14C W? ( 1 )=S( I3C ) 

W 3 ( 2 ) = S ( 13 ) 

CC TP 20 C 
15C W 3 ( 1 ) = C . 

U ? ( 2 ) = 1 . 

GP TC 2fC 
C 

169 I ? CC = I 3C 4-N T 

IF I I 3CC ) 19 C » 180, 1 7C 
C 

C 13 IN SFCOND QUADRANT 
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KB I T = K B I T/ ^ 

KL=KF IT-? 

C 

c nc tcp j=c 

cc ar i=i,iii,idif 

K L AS T = I + KL 
DC 3C K=I,KLAST,2 
K 1 =K +K B I T 
K2 = K 1 + KB I T 
K ?=K 2 + KBI T 

DC TWO STEPS WITH J = f 
A ( K ) =A (K ) + A ( K2 > 

A(K2)=A(K)-A(K2) 

A ( K 1 ) = A ( K 1 ) + A < K3 ) 

A ( K ? ) = A ( K 1 ) - A ( K 3 ) 

A(K)=A(K)+A(K1 ) 

A ( K 1 ) = A ( K ) -A ( K 1 ) 

A ( K2 )= A ( K2 ) + A ( K 3 ) * I 
A( K3) = A( K2 ) - A < K3 ) *1 

T = A ( K ? ) 

A ( K? ) = A ( K ) -T 
A ( K ) = A ( K ) + T 
T= A ( K? +1 ) 

A ( K 2 + 1 ) = A ( K+l ) -T 
A ( K + I ) =A ( K + l ) +T 

T = A ( K3 ) 

A(K3) = A(K1 )-T 
A ( K 1 ) = A ( K 1 ) + T 
T=A( K3+1 ) 

A(K3+1)=A(K1+1)-T 
A ( Kl + 1 ) = A( Kl + 1 ) +T 

T= A ( K 1 ) 

A ( K1 ) = A ( K ) -T 
A ( K ) =A ( K ) + T 
T=A( Kl+1 ) 

A ( K 1 + 1 ) = A ( K+l ) - T 
A(K+1)=A(K+1)+T 

R=-A( K3+1 ) 

T = A ( K 3 ) 

A( K3 )= A( K2 )-P 
A(K2)=A(K2 )+R 
A ( K 3 + 1 ) =A( K 2+ 1 ) - T 
80 A(K2 + 1 )=A( K2 + 1 ) + T 

IT (JLAST) 235,238,82 
82 JJ-JJPIF +1 

DC F CP J = 1 
I L AS T= II +JJ 
DC 65 I = JJ , HAST , ID I F 
KLAST = KL+I 
DO 85 K=I,KLAST,2 
K 1 = K + KBIT 

K? = Kl+KRIT 
K 3 = K 2 + KB IT 

LF TT ING W= ( 1 + I ) /RCCT2 ,W3= ( -1+ I ) /R0CT2 * W2 = I , 
A ( K ) = A ( K ) + A ( K2 ) * I 
A( K2 I =A(K) -A( K2 ) * I 
A ( K 1 ) * A ( K 1 )*W + A( K3) *W3 
A ( K 3 ) = A ( K 1 J*W- A( K3 ) *W 3 

A ( K ) = A ( K ) + A ( K 1 ) 

A ( K 1 ) = A ( K ) - A ( K 1 ) 

A ( K2 ) = A ( K2 ) + A ( K 3 ) * I 
A ( K3 ) = A( K2 ) - A ( K3 ) * I 
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XZ ERUXTPY 

IF (XZt RO.LT.0.0C 0 ^ 1 ) GO TO 5 
1 CONTINUE 
5 RETURN 

end 



SCPFC'UTINE HARM( A, Y, INV, S, IFSFT, IFERR) 

01 NF NS I ON A ( 1 ) , INVCl ) , S( 1 ) ,N( 3) ,M( 3) , NP( 0 ) ,W( 2) , W2( 2) , 
ECU VALENCE (N1,N( I)) ,(N2,N(2)),(N3,N(3)) 

10 I F ( IABS(IFSET) - I) 900,900,12 

12 Y T T= YA XC (Y(1),Y(2)»Y(3)) -2 
PC OT 2 = SORT ( 2 . ) 

IE (YTT-YT ) 14,14,, 13 

13 I T F F F = 1 
RETURN 

14 IFERR=r 
mi = R( 1 ) 

F 2 = M ( 2 ) 

N?=M(3 ) 

n i = 2 * * y i 

\ 2 =?**Y2 
N 3=2** M3 

16 IE(IFSET) 18,18,20 

18 NX = N1*N2*N3 
FI = NX 

cr i° i -- i , nx 

AC2M-1) = M2*I-i)/FN 

19 AC 2*1 > = -A( 2*1 ) / FN 

20 NP ( 1 ) = N1 *2 

N P ( 2 ) = N P ( 1 ) *N2 
NP(3)=NP( 2 )*N3 
CO 2 FT I D= 1 , 3 
IL = N P ( 3 ) -NP ( ID) 

I L 1 - I L + 1 

y i = y < i r ) 

IE (MI )25C ,253,3" 

3 0 IT IF =N PC ID) 

KP IT = N P ( ID ) 

MF V = 2 -*- (y 1 / 2 > 

IF ( M I - MEV ) 6C , 6C , 4 f 
C 

C v IS ODD. DO L =1 CASE 

40 KPIT=KBIT/2 
KC =KBI T-2 

DC 5C I=I,IL1,IDIF 
KL AS1 = KL* I 
DC 60 K= I , KL AS T , 2 
KT =K+K PI T 

DO ONE STEP WITH L=1,J=" 

A ( K ) = A ( K ) + A ( KD ) 

A ( KP ) = A ( K ) -A ( K D ) 

T = A ( KL* ) 

A ( KD ) = A ( K ) -T 
A ( K ) = A ( K ) + T 
T ■= A < KD 1 > 

AC KD+1 ) = A ( K+ 1 >-T 
50 AC K* 1) =A(K+1 ) +T 

IF ( MI - 1)250,250,5? 

52 L F I P ST =3 

C OFF - JLAST = 2**(L-2) -1 

J L AS T= 1 

GO TC TC 
C 

C y IS EVEN 

60 LF IRST = 2 
JL AST= 0 

7 0 DC 2 41 L =L E [ RST ,M I , 2 
JJCIE=KDIT 
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TAPF CONVERSION PROGRAM 



TH C ASSEMBLE R LANGUAGE PROGRAM 'FORM* THAT DOES THE 
ACTUAL PIT C ON VF F S l ON WAS WRITTEN BV P I MC.0 R N C.ZELENV 
CF ThF NAVAL POSTGRADUATE SCHOOL COMPUTER FACILITY STAFF 



C THIS PROGRAM CONVERTS 7 TRACK, 24 BIT WORD TAPE 
C TO S TRACK 22 SIT WORD TAPF. 

DIMENSION INDATAI 1024) , DAT A ( l '24 ) 

F ACT PR = Iff. . V( 2**U-1 ) 

REWIND 2 
REWIND 4 
NC BL K S = 1 

C DO J= 1 , K WHERE K IS THE NUMBER OF BLOCKS OF 1024 EAC 
Du 31 J=1,I CO 
READ (2,3, FRR=5C ) INDATA 
CALL FORM! INDATA) 

DO 1 1=1,1024 

1 DATA ( I ) = INiTATA ( I ) -‘FACTOR 
WPITF (4,^) DATA 
CALL P MEAN! ( 0 AT A » N OBLK S , 8MI » PM S ) 

W R I T F ( 6 » 7 0 )J,BM,QMS 
GL TO 31 

d "* WPITF (6,51) J 
31 CONTINUE 
3 FORMAT ( 1 6 ( 64A4) ) 

51 FI RMAT ( ,r, ,5X,'READ ERROR, RECORD N0. = ',I6) 
ofc FORMAT ( IX, 1 ^ E 1 C . 2 ) 

72 FORMAT ( 1CX RECORD NQ=',I6,' MF. AN = * , Frt . 3 » ' MEAN SQ 

* F 1 0 . 'a ) 

STOP 
Of D 



SI "ROUTINE 3 M F A N ( P , NO HL KS , BM , BMS ) 
DIMENSION P(l) 

MF =0 
B M S = ' .0 
BM = ”' .C 

1 D( . 5 1=1 ,NiCBLKS 
J=( 1-1 KT 24 
SOM=f .R 

DC 3 I 1 = 1 , 1 '24 

J J =J + I I 

IF ( M K . FO. 0) GO TO 2 
SLM=SUM+FM JJ ) * *2 
G(' Tf' 3 

2 SOM=SUM+P( JJ ) 

3 CONTINUE 
SGM= SUM / K 2-t- 

II (MK . EG. ;• ) GO TO 4 

BM S=P>MS+ SUM 
GO TO 5 

4 HM=P.M+SUM 

5 CONT INUO 

I F ( MK .EC. ^ ) GO TO 6 

PMS=BMS/NOPLKS 

Gi TO fi 

6 P M = P M / NO P L K S 

DC 7 1=1, JJ 

7 P ( I ) = P ( I ) - PM 
MK =1 

GO TO 1 
3 CENT IN LIE 
RETURN 
ENC 
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SliBFt’liT Ih F F ( i R M ( INCATA ) 

F f] •< M FTAPT ' 

* SUIRCUTINF FCP41INDATA) 

* 

* THIS SUBROUTINE SILL CCNVFRT 24 BIT BINARY WOKCS STORED IN 

* ^N ARRAY LFROTF SPECIFIED 6 Y THE INDEX VALUE TO M HIT PIN 

* U.C PLACE THESE SAME WORDS HACK INTO INDATA 



* 

* 



LOOP 



CAT A 
NUM 



STM 14,12,12113) This subroutine converts 



E/LR 


6, j 


24 


HIT 


BINARY WORDS TP 


US INC 


* , 6 


22 


B I T 


WORD S 


U S I KG 
SR 

l 


DATA, 7 
7, * 

1 1,=F »1C 24* 


THIS 


I S 


THE INPEX VALUE 



L 1 2 , r ( 1 ) 

L 2,NUM(12) 

L u 2,7 
SR PL 2,6 
SEL 2,2 
SE OL 2 , j 
SR L 2,2 
SR PL 2 , *> 

SFL 2,2 
SPEL 2,6 
ST 3 , NUM (12) 

LA 12,4(12) 

ECT 11, LOOP 

L M 2 , 1 2 , 2 0 I 1 3 ) 

MV I 1 21 1 3 ) ,X » FF* 

BCR 15,14 

DS ECT 

DS 1 F 

FT P 
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DISPLAY PROGRAM 



TH 11 SUBROUTINES THAT GIVE THE POINT PLOT ON THE I3M 2 25C 
DISPLAY unit WERE WRITTEN BY LT.(JG) G.J. VORHUFF, 
DEPARTMENT (F AF RON AUT I CAl ENGINEERING* 

NAVAL PL STOP APUA T E SCHOOL 



INTEGER HOLE 

DIME NS ION Y< 5~* ) ,CRU(2b) *CPL( 25 ) , RANGE ( A) , ETA (25 ) * 
*CEM<?,?5> 

R F A L * a TITLE! IS)/' UNSTEADY PRESSURE COEFFICIENTS VS ET 
* ' / 

NO I M= 2 3 
RANGE' ( 1 ) = 1 . 

RANGE ( 2)=C . 

RANGE ( 3 ) = 6 . 

RANGE <4)=-l. 

KE AO( 5, If I ) ( ETA( I ) , 1 = 1 ,25) 
i : 1 FERMAT ( 16F5.C ) 

CALL GOGSP 
9 597 CL N'T INUE 

DC 1 I T= 1 , 1 24 
RE AD <4,1: 2) Y 
K'2 Ft RM AT ( 5C A 4 ) 

IF ( IT .NE. I ) GO TO 3 
l =C 

DO 2 1=2, .SC, 2 
L = L+ 1 

CEL<L)=-Y( I - 1 ) 

C P M ( 1, L ) =C PtJ ( L ) 

CPU L) = -Y( I ) 

2 C PM ( 2 , L ) =C PL < L ) 

GL TO 6 

3 DC A I =1 ,2 5 

A CPU( I ) =-Y< I ) +C P w ( 1,1) 

DC 5 1=26,5 
J= 1-25 

5 CPU J)=-Y( I ) +C PM ( 2 , J ) 

6 C( NT INUF 

CALL LI TPC R T ( ET A, CPU » NO IM , RANGF , 1,1, TITLE, £9998, £99 99) 
CALI UPHC R T ( ETA,U»L,NDI GRANGE , 1,3, TITLE, I J, I K , £ 9 9 9 8 , £ 
1 CONTINUE 
9*)8 CONTINUE 
CALL CGSP 
GC TO 9997 
9959 CONTINUE 
STOP 
EE C 



SUE PC UTI NE UTPCRT ( U , V , NOATA , RANGE ,KKZ , MOOCUR , TITLE,*,* 

UTPCET PrOCUCFS A PUINT-PLOT ON THE 225C GRAPHIC DISPLAY. 
U IS THE VECTOR OF ABSCISSAE 
V IS the vrCTCU OF ORDINATES 

NOATA IS THE LENGTH CF VECTORS 0, AND V (THE NUMBER OF P 
RANGE IS THE SCALING VFCTUR: 

RAECE(l) = MAXIMUM DESIRED X VALUE PLOTTED 

RANGF (2) = MINIMUM DESIRED X VALUE PLOTTED 

RANGE ( 3 ) = MAXIMUM DESIRED Y VALUE PLOTTED 

RAEGLIM = MINIMUM DESIRTD Y VALUE PLOTTED 

EVERY KKZ'TH ELEMENT OF THE INPUT LISTS WILL BE PLOTTED 
USER SETS MCDCOR = C 

TITLE IS EITHER THE NA^E OF A 6 ? - BY T E ARRAY UR A 6'.-CHARA 
CONST ANT. 

DATA SET FTC9F'M MUST BE PROVIDED, DC 8= ( RFC F M = F , LR E C L= 8C 
FIPST RETURN IS FOP RESTART, SECOND FOR TERMINATION 



C 



COMMON/ IGNORE/ NULL , I DEV, IGSPNM 
D I MR NS IOE A ( 4 ) ,6(4) , XSCALF( 5) , YSCALE( 5) 
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1 t MJL L ( 1 ) ,APEA( 19 ) , T ITLEU5I 
o i me ns i un u< i) ,v( i ) ,x( ia ) , y< nr ) ,range(M 

RE AL* 4 XCHAP ( 4 ) / 1H+ , 1H* , IHX , 1H. / 

1 1' TEGE k A, B 

if (Monr.uF . eg. ' )wute(6,4| title 

4 FT PM AT < •< • , 1 5A4/// ) 

CALL U TP LPT ( U , V , N D AT A , R ANGE , KK Z , MOOCUR ) 

DP 2 I = 1 » f i DA T A 
X ( I )=U< I ) 

Y ( I ) = V ( I ) 

2 CPNTINUF 

INITIALIZE GRAPHIC DATA SETS 



CALL I NODS ( IDEV, IGOSNM) 

CALL S DAT M ( IGOSNM , 1 ) 

CALL SGDSL ( IGOSNM, l.P , 1.64, 1 1 .3 , 1 1 .5 , r . , 0 . , 1 2 . r L 2 . ) 

X M AX = RANGE ( 1 ) 

X M I N= RANGE (2 ) 

YMAX=P ANGE ( 3 ) 

VM I N = R ANGF ( A ) 

XR ANt.E = XMA X-XM I N 
YE Af GF = YMAX-vmIN 
DX=XRA\CF/ 53. 
f)Y=YRANGE/M . 

X L L = XM IN-DX 
XUR=XMAX+DX 
YL L= YM I N-D Y 
VU«=YMAX+DY 

CALL SDATL ( I GDSNM , XLL , YLL , XUR, YUR ) 

CALL I N(,DS( IDEV, IGDSNN) 

CALL SLATM (IGDSNN ,3) 

CALL SDATL (IGDSNN ,0,',73,Sl) 

PLCT TjTLE 

CALL PTEXT (IGDSNN , TITLE, 6 ‘ , NULL » NU L L » NU LL n » 5 l ) 
CALL F X F 0 ( IGDSNN, £9998, €9999) 

PLCT GRAPH T R AME 



A ( 1 ) = 1 1 
B( 1) =7 
A ( 2 I = 1 1 
6(2)=^° 

A ( 3 ) =t 9 
H( 3 ) = 49 
A ( 4 ) =o 9 
B( 4) =7 

CALL STPOS ( IGDSNN ,69,7) 

CALL PLINE( IGDSNN , A , B , NULI , MULL , NUL L , 4 ) 

CALL EXFC( IGDSNN, £9998, £9999) 

CCMPUTE PR( PEP SCALE NUMBERS 

3 j C 0 X INCR=XPANGE/4. 

Y I NCR’ = YP ANGE/4 . 

XSCALE ( 1 )=XMAX 
YSCALE ( 1 )= YMAX 

or ec 1=2,5 

9? X SCALE ( I ) = XSCALE ( 1-1 1 -XI NCR 
DC 8 1 1 = 2, 5 

81 YSCALE( I ) = YSCALF( 1-1 I-YIMCR 

PLlTTING X AND Y AXIS, IE NECESSARY 

YTEST=YMAX * YMI N 
XT EST=XMAX «XMI N 
If (XTEST) 1 ,222,222 
222 IE (YTEST ) 333,444,444 
1 CPNT I NUE 

CALI STPPSUGDSNM ,E.,YMIN) 
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CALL 


PUNE ( IGOSNM 


, a .,YMAX) 


CALL 


EXFGI IGOSNM, 


£9«98 , £999 ^ ) 


GOTO 


222 




333 CCNTINUF 




CALL 


STPOS ( IGDSNM 


, X Y I N , 0 . I 


CALL 


PL INF ( IGDSNM 


, XMAX, C. ) 


CALL 


FXFGI IGDSNM, 


£9998, £9999 ) 


444 CONTI NUF 





c 

c pli:t scale lumbers 
c 



REWIND 9 

DC 1 C 2 1 = 1 ,5 

WRITE ( c » 1 P ) V SC ALE ( I ) 

18 FORMAT < G 1 ' . 3, • - • ) 

10? CONTINUE 

wPITF(9,5M 

64 FERMAT ( 1 2X , • | • t A ( 1 3 X , • | • ) ) 

WRITE ( ° , 2 2 1XSCALEI 5) » XSCALE(4),XSC4LE<3),XSCALF(2),XS 
22 FORMAT ( 4 X , 5 ( 4 X » G 1 ".31 ) 

PNC TILE 9 
REWIND 9 
CP 181 1=1,5 

READ! 0 , 53) (AREA(K)»K=1,19) 

33 F( PM£ T ( 1 c A 4 ) 

L INE=43-( I - 1 ) *i: 

CALL PTEXT ( IGOSNN , AR E A , 12, NULL , NUL L , NUL L » 0 » L I NE ) 

Cm LI E XFC( IGDSNN, C999B ,£99QQ) 

1'1 CENT I N U f 

DC: 1 C 3 1 = 1,2 

P. FAD(°,53) (AREA(K) ,K=1,19) 

L I NF = P - 1 

CALL PTEXT ( I GDSNN , A REA, 74, NULL, NULL, NULL, 3, LINE) 

CALL E X E 0 ( IGOSNN, £9993, £9999) 

1 13 CONTINUE 

I TEMPI =IGDSNM 
I T EM P2 = I GD SNN 
If RR=C 
GC TO 5 
C 

C (NOTE: FNTPY UPHCRT ALLOWS FOR ADDING CURVES TO GRAPH WH 

C HAS AIFEADY BEEN LAID OUT -- CALLFD BY PRNCRT) 

C 

FNTPY UPHCPTUI,V,NCATA, RANGE, KKZ, MCDCUR, TITLE, IGDSNK, I 
C 

ir (MCOCUK. EM. 3 )WR l TEC 6,4) TITLE 

CALL UTHOLDI U, V, NDAT A , RANGE , KKZ , MODCUR ) 

D(J 3 I = 1 , N D A T A 

x ( n=u ( 1 1 

Y ( I ) =V ( I ) 

3 CCNTINUF 
5 CCNTINUF 

I GCSNM=ITF MP1 
I LDSNN'= I TE MP2 
C 

C CHECKING X AMD Y POINTS AND PLOTTING THOSE OUT OF RANGE A 
C 

DC ?l 1=1, N DATA, KKZ 
I E ( X (II-XMAX) 2 15,?'. 5,220 
??8 X ( I ) =XMAX + DX 
I F PE = I FRP + 1 
GC TO 2 1C 

2'-5 I E ( X U)-XMIN)2' 3,210, 211 
?.j 3 X ( I ) =X M I N -DX 
IFPP= IH-R + 1 

213 If (Y ( I J-YMAX 1215,215,212 
212 Y ( I > = YMA X ♦ D Y 
I F R R = IFFR + 1 
GC TO y 

>15 IF (Y (I ) - YM I N ) 2 1 7 , 33,33 

2 1 7 Y (I ) =YMI N -DY 
IFRF = IfRKM 
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30 CCNTINUF 
C 

C PLACING POINTS IN THEIR PROPER POSITIONS 
C 

IF <M( DOUR. F J . r ' . Ck . MC5DCUR . FO. 1 ) JSFT = ' V 

J S FT = J SE T + 1 

If ( JSET.GT .A) J S E T = 1 

DO 91 1=1 , NUAT A , KKZ 

RXCHA*=XCHAR( JSFT ) 

X I =X II ) 

V I = V < I ) 

CALL PTEXT ( I GDSNM , RXCHAR , 1 , NUL L , NUL L , NUL L , XI , YI ) 
CALL F X L 0 ( IGDSNK,£999«, 69999) 

CE NT I Ml F 

I F ( PC.P CUR . t J. " . C° . MOOCUR . FO . 3 ) GO TO 9?.?. 

PE TUP F 

T22 CONTINUE 
C 

IF ( I F P F. ) 1 C V>, 1" * " , 1 1 

10 0 1 CENT I NUL 
REWIND c 

WR ITE( 9,2c ) IPRR 

2~ FERMAT ( 'NUMBER OF POINTS OUT OF RANGE = ' , I A ) 
ENT FILE ° 

REWIND 9 

R I AD ( a , 5 3 ) (APEA(K),K=1,19) 

CALL PTEXT ( I GUSNN , AREA, A , NULL , NUL L , NUL L , ° » A ) 
CALL F XCU( IGDSNN|C9998,C9999) 

1 ):* CONTINUE 

C IE (MLTCDF . EG. 3 IKETUPN 

C 

C WAIT FUR END-KEY ATTENTION 
C 

CALL C P E A T N ( I A T T N 1 

CALL PC ATT ( I ATTN, I CODE , NULL , 699 93, 69999) 

CALL NPATL ( I ATTN, -I ) 

CALL F NATL ( I AT TN ) 

CALL T MGCS ( I GDSNM ) 

CALL TMGDSI IGDSNN) 

PC TURN 

9993 CLNTIFUF 
R f TUP N 1 
9999 CINTINUE 
RFTUFN 2 
END 



SUBROUTINE PRNCPT ( X , Y , NN , MODC UR , U TL E , * , * ) 

C 

C PRNCPT PRODUCES A PCINT-PLGT UN THE 2 25'' DISPLAY. IT CAN 
C SEVFFAI CURVES ON THE SAME GRAPH, DEPENDING ON THE CALLED 
C MCDCGR: 

C , Tf-IS IS THE ONLY CURVE FOR GRAPH 

C =1, THIS IS THE riRST jf SEVERAL CURVES FOR GRAPH 

C =2, THIS IS AN INTERMEDIATE CURVE FOR GRAPH 

C =3, THIS IS THE LAST OF MANY CURVES FCR GRAPH 

C X IS THE VICTOR OF ABSCISSAE 
C Y IS THE VFCTOR OF ORDINATFS 

C NN IS THE LENGTH OF VECTORS X AND Y (THE NUMBER OF POINTS 
C TirLf IS EITHER Thf NAME OF A S^-BYTE ARRAY OR A 6P-CHARA 
C CONSTANT. 

C PRNCPT CREATES THE PROPER RANGE VFCTOR AND CALLS UTPCRT 
C PCIMS OF THE FIRST CURVE ON A GRAPH. SUCCEEDING CURVES 
C ARE TRUNCATED TO THF SA^T RANGE AND PLOTTED WITH DIFFEREN 
C FIRST RFTUPN IS FOR RESTART, SECOND FOR TERMINATION 
C DATA SET EK9F '1 RUST EE PROVIDED, DC B = ( RE CF M= F , LR F CL = 8 2 
C 

CC MM (It. / 1 GNPq F / NULL » I DEV, IGSPNM 
C 

DIMENSION X(t 1), Y(5I), RANGE(A) 
l ULL ( 1 ) , APE A ( 19 ) , T I r L E ( I 5 ) 

ECUIVALENCF ( R ANG P ( 1 ) , X'1 AX ) , ( R ANGE ( 2 ) , X M I N ) , ( RANGE (3) 
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i 



( R A N G E ( 4 ) | YM I N ) 



BRANCH CN MODC'JR 

IF (MCDCUR . EG.' 1 . CR . MQDCUR . GQ . 1 ) GO TO 40 
GL TC 4'-* 

G AUG T INPUT DATA AND COMPUTE RANGE FOR FIRST OR ONLY CURV 
4? XNAX=-l.F2r 



X F I N= 1 . R2 r 
YMAX--1.P2C 
YMIN=1 .E2C 

nr i i = i » n n 




1 F (X ( I l-XMAX ) 
XMA X = X ( I ) 

Y X m A X = Y ( I ) 


6 , 6 , 2 


IF (X ( I 1-XMIM) 
XM I N=X ( I ) 

YX MIN=Y( I) 


3,3,7 


IF (V( I ) -YM AX ) 
Y M A X = Y ( I ) 

XV MAX=X ( I ) 


3,8,4 


I F ( Y ( I ) - YM I N ) 
YM IN = Y( I ) 

XV M I N=X< I ) 

Cf NT IF OF 


5,5,1 



CALL LTPCKT PROPERLY 



4 r . r IT (RTtCUR. EQ.2.CR.MC0CUR.EQ.0) GO TO 

CALL UTPCRTI X , Y ,NN .RANGE , 1 , MUOCUR , T I T L E , L 9 9 9 8 , C999 9 ) 

GL TO 7C 

50 CALL UPHCRTI X, Y,NN, RANGE, 1 * MCDCUR » T I TL E » I G C , I G01 , 6999 8 
70 I r- (MCDCUR .NE . 3 IRFTURN 



LIST TRUNCATION RANGE AFTER END OF MULTI-CURVE PLOT 



' 1 



61 



6 r - 



REMIND 

V. R I T F ( 9 , 1 ' 1 ) Y PA X f 
FORMAT ( 'MAX Y = • 

1 • M N Y = » 

NR ITE ( 9, im X M A X » 
FORMAT ( 'MAX X = ' 

1 'MIN X = • 

ENC RILE R 
I C CR NM = I Gf 1 
REFIN'D 9 
DO 60 I = ac t 52 
RE AD ( c ,6 I ) 

FORMAT ( 1 c A 4 ) 

J= 92- I 



XYMAX t YMI Nt XYMI N 
, G 1 4 . 7 , • AT X = • ,G 1 4. 7/ 
, G 1 4 . 7 , ' AT X = • , G 1 4 . 7 ) 
YXMAX , XMIN, YXRINi 
, G 1 4 . 7 , ' AT Y = ' , G 1 4 . 7 / 
,G14 .7, • AT Y = • ,G14. 7 ) 



(APEA(K) »K=lt 19) 



CALL PTLXTIIGDSNM, AREA, 74, NULL, NULL , 1 ,0, J ) 
CALL R XRQ( IGOS N'M * 6 99 Ci 8 , 69999 , 0.62 ) 

CONT I NOE 



WAIT FOR END-KEY ATTENTION 
CALL CREATM IATTN ) 

CALL RCATTI I ATTN, I COO E , NULL , 69998 , 69999 ) 
CALL MPATL ( IATTN, -1 ) 

CALI RNATL(IATTN) 

62 CONTINUE 

CALL TMGOSMGO) 

CALL TMGDS(IGOl) 

Rf TURN 
53 CONTINUE 
RETURN 1 
9979 CONTINUE 
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RETURN 2 
END 



SUBROUTINE i)ELAY(N) 

CALSF TINT LAPSE OF N HUNDREDTHS OF A SECOND (REAL TIME) 

INTEGER T1PE,0T 
T I Mt = I TIKE (0 ) 

1 CONTINUE 

0 T = I T I P t (D)-Tl M E 
IF(DT.LT.N ISO TO l 
RE TURN 
END 



SI PRC'U TINE tXEQ( I GCSNM, *,*,*> 

C OF T R CL DISPLAY REGENERATION — HOLD FURTHER GRAPHIC DATA 
BUTTER IS P F A C v 

COMMON// TTFNT/IATNL 

C C PPCN/ I GNLRE / NUL L » IDEV, IGSPNM 

CALL I- CATN( I ATNL , I CODE , I , NULL ,0 , 3 1 ,32 ) 

IF ( ICCDF.F C.-l ICALL START ( &R^3) 

I r II CODr . t C. 31 ICALL F I N I S H ( £^99 8 , 09OQ9 ) 

If < ICLDE .FG. 32 I PE TURN 3 
CALL FXEC(IGnSNM) 

CALL f f L AY { 5 ) 

RETURN 

' Cl. NT I NUF 
RFTUFN 1 
3-J9 CO IT I NUF 
RFTURM 2 
ENT 



SUBROUTINE GOG SR 

INITIALIZE CSP, 2 25 r , AND SEVERAL WORK-DATA-SETS 

CfPPCN/ATTtNT/ I AT NL 

CC OMON/ IGNORE /NULL »inPV,IGSPN M 

N U LL = - 5 

CALL I NGSP ( IGSPNM, NULL ) 

GC TO. 1 

cNTRY PROVIDES FOR P E - I N I T I AL I L AT I ON 

ENTRY GGSP 
CALL T PDF V ( IDEV ) 

l C Of T I NUF 

CALL INDFV ( IGSPNM, 22» IDEV) 

CALL S A L f- P ( IDEV) 

CALL PL I TS ( I HE V, 3 } 

CALL f RATL ( I DL V, I ATNL ) 

CALL ENATNd ATNL, 3, 31 ,32) 

RETURN 

en r 



SLOP OUT INF START!*) 

RESPLNSF TO R E- I N I T I AL I Z A T ION COMMAND DURING EXECUTION 

C CPPON / IGNORE/ NULL , I DEV , I GSPNM 
C 

C A l L SALPM (IDE V ) 
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CALL TNCFVL I0EV) 

CALL I NOEVI IGSPNM,22, IDEV) 

CALL INITUSI IGCPGL,3,2,' 1 .,0. , 1 . , 1 . , ^ . , ? . , 4 8. , 34 . , . TRUE 
DC 1 1=1, 1C 

CALL PTEXT ( IGDPGL , 'RESTARTING G SP . . . ' , 17 , NULL , NUL L , NUL 
1 CONTINUE 

CALL EXEC ( IGDPGL ) 

CALL D EL A V ( 1 •**'? ) 

RETURN 1 
E N 0 



SUBROUTINE. FINISH!*,*) 

RESPCNSE TO TERMINATION PECUfc'ST DURING EXECUTION 
( ALL l A; ONE CHANCE TO CHANGE MIND) 

COMMON/ IGNORE /NOLL , I DEV, IGSPNM 

DIMENSION FMTI13) 

LATA FMT/ • ( • • 1********** EXECUTION ABORTED 3Y USER *** 
CALL SALRM ( IOEV) 

CALL TMDEV(IDEV) 

CALL INDFV( IGSPNM, 22, IDEV) 

CALL I MTDSI IGDPGL,^,?, , 0 . , 1 . , 1 . , 3 . , 0 . , A 8 . , 34 . , . TRUE 

INTENSIFY MESSAGE 

CALL 8GSE Q (IGOPGL ) 

DC 2 1=1, 1C 

CALL PTEXT ( IGDPGL ,• LAST CHANCE... START OR FINISH?', 31 
1NOLI *9,17) 

2 CONTINUE 

CALL ENSEWI IGDPGL , KEY ) 

CALL E XEC ( IGDPGL ) 

CALL CRATL ( I DEV, I ATNL ) 

CALL F N A T N ( IATNL, ',31 ) 

F L AS F FNKB LIGHTS AND CRT DISPLAY 

1 CCNTIMJE 

CALL DELAY (5^) 

CALL CMI T( IGDPGL , NULL ,KPY) 

CALL ML I T S ( IATNL, 2) 

CALL DEI AY (GC ) 

CALL I NCL ( IGOPGL, NULL , KEY) 

CALL ML I TS ( I ATNL , 3 ) 

CAIL PCATNI I ATNL , ICODF,l,NULL ,0,-33) 

I F ( I CODE. F G. 3) GO TO 1 
CALL ML ITS ( IATNL , 1 ) 

I F ( JCC0E.FC.~1 ) C AL L START (£99Q8> 

CALL R CSFT ( IGDPGL ) 

C CNF IRM TERMI NATION 

CO 3 1=1, 1C 

CALL PTEXTUGDPGL ,FMT(2) , 44, NULL, NULL, NULL, 3,17) 

3 CONTINUE 
CALL E X EC ( IGOPGL ) 

V,FITF(6,FMT) 

CALL DELAY (G^D) 

CALL TMDEV(IDEV) 

RETURN 2 

9993 CONTINUE 
PFTUPN 1 
CND 
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noon o o ooo 



SUBROUTINE INI TDSI IGDSNM, 1DATM, l CH AM , X l , Y l , X 2 , Y ?, X3,Y3 
I NT F F F A C L TC GSP GDS -CRE AT ION AND OPT l JN-OEF l NI T I ON ROUTI 
CC MM ON / I GN ORE / MLL L • IOEV, IGSPNM 
LOGICAL LP 

CALL INGDSI IDEV, IGDSNM) 

CALL SGI' SI. ( IGDSNM, XI, Y1,X2, Y2 ,G. , , 1 . ,1. ) 

CALL SDATLI IGDSNM, X3,Y3,X4,Y4) 

CALL SDATMl IGDSNM, IDATM) 

CALL SCHAM.I IGDSNM, ICHAM) 

IF ( LP ) CALL SLPATI IGDSNM, 1) 

RF TURF 
El C 



SUBROUTINE C RE ATN ( I AT TN ) 

INTERFACE TC GSP 'CRATL' 

MAINTAIN RFSTAP T/ABORT LEVEL AS ACTIVE 

Cf R vr JN/IGNCRE/NULL » I DEV, IGSPV* 

CALI C RATL ( IOEV, I ATTN ) 

C A I L MLPFC(IATTN,1,2,1) 

CALL NLPEOI I AT TN , 2 , 4 , ] ) 

CALL MPATL ( I ATTN, 1) 

CALL FNATM I ATTN,0 ,31 ,32) 

PF TURN 
END 



SUBROUTINE PQATT ( I ATTN, I CODE , / I RR AY / , «■ ,* ) 
INTERFACE TC GSP •PQATN' 

RESTART OP ABORT ATTENTIONS TAKE PRECEDENCE 

CCMMGN/A1 TENT/ I A TNL 
COMMON/ IGNORF/MULL , IDEV, IGSPNM 
C 

DIMENSION IRRAY(l) 

CALL PCATN ( I ATNL , I CODE, I , NULL ,0 , 31 , 32 ) 

I F ( I CODE .EC. -I ) C A L L START ( £9990 ) 

I F ( ICODF.FQ. 31 )CALL F I N I SH ( £9998 , £9999 ) 

IF ( I CODE . F 0 . 32 ) RETURN 
CALI MPATL ( I AT TN » -1 ) 

CALL R GA TN ( I AT TN , I CODE ,2 , I RRA Y , 9 , - 35 ) 

CALL MPATL U ATNL , -1 ) 

IF ( ICCJDE.CC.-l )CALL START (£9998) 

IF ( ICC0E.F0.31 ICALL F l N l SH ( £999 8 , £ 9 999 ) 

PE TURN 

999R CONTINUE 
RETURN 1 
9999 CONTINUE 
RETURN 2 
END 
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